TRANSONIC SHOCKS IN MULTIDIMENSIONAL DIVERGENT NOZZLES 

MYOUNGJEAN BAE AND MIKHAIL FELDMAN 

Abstract. We establish existence, uniqueness and stability of transonic shocks for steady 
compressible non-isentropic potential flow system in a multidimensional divergent nozzle with 
an arbitrary smooth cross-section, for a prescribed exit pressure. The proof is based on 
solving a free boundary problem for a system of partial differential equations consisting of 
an elliptic equation and a transport equation. In the process, we obtain unique solvability 
for a class of transport equations with velocity fields of weak regularity (non-Lipschitz), an 
infinite dimensional weak implicit mapping theorem which does not require continuous Frechet 
differentiability, and regularity theory for a class of elliptic partial differential equations with 
discontinuous oblique boundary conditions. 



1. Introduction 
We consider inviscid compressible steady flow of ideal polytropic gas. 

One of important but difficult subjects in the study of transonic flow is to understand a global 
feature of flow through a convergent-divergent nozzle so called a de Laval nozzle. According to 
the quasi-linear approximation in |CFR1 Chapter V. section 147], given incoming subsonic flow 
at the entrance, the flow is always accelerated through the convergent part of the nozzle unless 
the exit pressure exceeds the pressure at the entrance. The flow pattern through the divergent 
part, however, varies depending on the exit pressure and the shape of the nozzle. In the 
divergent part, the flow may remain subsonic all the way to the exit, or it may be accelerated 
to a supersonic state after passing the throat of the nozzle and have a transonic shock across 
which the velocity jumps down from supersonic to subsonic. In particular, the approximation 
implies that if the exit pressure is given by a constant pc satisfying < Pmin < Pc < Pmax < oo 
for some constants Pmin and Pmax^ then a transonic shock occurs in the divergent part of the 
nozzle. One may refer [CFR] for the quasi-linear approximation in curved nozzles. In section 
2.41 for any given constant exit pressure G (PmimPmax), we rigorously compute the transonic 
flow of the Euler system in a multidimensional nozzle expanding cone-shaped. 

The issues above motivated recent works by several authors on existence and stability of 
steady transonic shocks in nozzles using the models of potential flow or compressible Euler 
system [CFTl [CF2l [CF3l IXYTl ICHTl lYUTl [CCFllXY2] . In these works transonic shocks were 
studied in cylindrical nozzles and its perturbations. A physically natural setup for transonic 
shock problem is to prescribe the parameters of flow on the nozzle entrance and the pressure on 
the nozzle exit. In particular, an important issue is to show that the shock location is uniquely 
determined by these parameters. However, such problem is not well posed in the cylindrical 
nozzles. Indeed, the flat shock between uniform states can be translated along the cylindrical 
nozzle, and this transform does not change the flow parameters on the nozzle entrance and 
the pressure on the nozzle exit. This provides an explicit example of non-uniqueness. Also, for 
uniform states in cylindrical nozzles, the flow on the nozzle entrance determines the pressure 
on the nozzle exit for a transonic shock solution. This degeneracy leads to non-existence, i.e. 
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one cannot prescribe an arbitrary (even if small) perturbation of pressure on the nozzle exit. 
On the other hand, as the de Laval nozzle example suggests, it is natural to study transonic 
shocks in the diverging nozzles. This removes the translation invariance in the case of Euler 
system. Recent works by Liu, Yuan [LY] . S.-X. Chen |CH3j . and by Li, Xin, Yin |LXYj showed 
well-posedness of transonic shock problem with prescribed flow on entrance and pressure on 
exit of a divergent nozzle, which are small perturbation of the corresponding parameters of the 
background solution, in the case of compressible Euler system in dimension two, under some 
additional restrictions on the perturbed pressure on the exit. In particular, it was assumed that 
for the prescribed pressure on exit, the normal derivatives on the nozzle walls vanish. While 
this condition is not physically natural, mathematically it leads to substantial simplifications: 
it allows to show regularity of the shock up to the nozzle walls, which allows to work with 
vector fields in the transport equations. 

In this paper we study transonic shocks in diverging nozzles in any dimension, for the 
case of perturbed cone-shaped nozzle of arbitrary cross-section. Moreover, we do not make 
assumptions on the pressure on the nozzle exit, other than the smallness of the perturbation. 
This requires to consider shocks with regularity deteriorating near the nozzle walls, and then 
vector fields in transport equations are of low regularity (non-Lipschitz). 

Furthermore, we study this problem in the framework of potential flow. A surprising feature 
of this problem is that it is not well-posed for the "standard" potential flow equation. Namely, 
for radial solutions in a cone-shaped nozzle, the spherical transonic shock can be placed in any 
location between the nozzle entrance and exit, without changing the parameters of flow on the 
nozzle entrance and exit. This follows from a calculation by G.-Q Chen, M. Feldman |CFlt 
pp. 488-489], see more details below. Thus, spherical shocks separating radial flows in the 
cone-shaped divergent nozzles can be translated along the nozzle for potential flow equation. 
This is similar to the case of flat shocks separating uniform flows in straight nozzles. 

In order to remove this degeneracy, we introduce a new potential flow model, which we call 
non-isentropic potential flow system. It is based on the full Euler system (as opposed to the 
"standard" potential flow equation, which is based on the isentropic Euler system), and allows 
entropy jumps across the shock. The system consists of an equation of second order, which 
is elliptic in the subsonic regions and hyperbolic in the supersonic regions, and two transport 
equations. 

The main purpose of this paper is to establish the existence, uniqueness and stability of a 
transonic shock for non-isentropic potential flow system in a multidimensional divergent nozzle 
with an arbitrary smooth cross-section provided that we smoothly perturb the cone-shaped 
divergent nozzle, incoming radial supersonic flow and the constant pressure at the exit. 

In Section [21 we introduce the non-isentropic potential flow model, and define a transonic 
shock solution of this new model. Then we state two main theorems of this paper and present 
a framework of the proof for the theorems. Section [3] is devoted to a free boundary problem 
for the velocity potential if where the exit normal velocity is fixed . In section [H we plug 
a solution ip of the free boundary problem in section [3] to a transport equation to find the 
corresponding pressure p. In section [Sj we prove the existence of a transonic shock solution of 
the non-isentropic potential flow for a fixed exit pressure pex- The uniqueness of the transonic 
shock solution is proven in section [6j 
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2. Main theorems 

2.1. Compressible Euler system and shocks. We study transonic shocks of steady inviscid 
compressible flow of ideal polytropic gas. Such flow is governed by the conservation of mass, 
momentum and the conservation of energy. These three conservation laws provide a system of 
PDEs for the density /), velocity u = (ui, • • • , Un) and the pressure p as follows: 

div{pu) = 

div{pu (S> u + pin) = ^2 1) 

div{pu{^\u\^+ ^^^_^^^^ )) =0 

where I„ is n x n identity matrix. The quantity c, given by 

(2.2) 




is called the sound speed. The flow type can be classified by the ratio, called the Mach number, 
of |u| to c. If > 1 then the flow is said supersonic, if < 1 then the flow is said subsonic. 

If \jA = I then the flow is said sonic. Depending on entrance and exit data of flow or a shape 
of a channel, there may be a jump transition of the flow type across a curve or a surface. Such 
a transition can be understood in a weak sense as follows. 
For an open set Q, C M", if {p,u,p) G L\^^{Q) satisfies 

/ pu-V^dx= I [pukU + pck) ■ V idx = I pu(^-\^'^ + ——^^^-—) ■ S/^dx = (2.3) 

for any ^ € Cg°(r2) and all k = 1, • • • ,n where is a unit normal in k-th direction for M", 
{p,u,p) is called a weak solution to ()2.ip in Vt. 

Suppose that a smooth (n-l)-dimensional surface S divides into two disjoint subsets il^ . 
If {p,u,p) G Ljoci^) is in C^{^^) C°{^}^) then {p,u,p) is a weak solution of ([2TT]) if and 
only if {p,u,p) satisfies (j2.ip in each of 0^ and the Rankine- Hug oniot jump conditions {which 
is abbreviated as R-H conditions hereafter) 

[pu-Ks]s = (2.4) 
[p{u ■ Us)u + pi^s\s = (2.5) 

for a unit normal Vg on S where [F]s is defined by 

[F{x)]s:=F{x)\-^-F{x)\-^ for x G 5. 

To study a transonic shock in M" for n > 2, we consider an irrotational flow in which the 
velocity u has an expression of = Vip for a scalar function (p. Such a flow is called a potential 
flow and tp is called a (velocity) potential. 



2.2. Steady isentropic potential flow equation and radial transonic shocks. A widely 
used steady potential flow model(see e.g. [CFHIXYI] et al.) consists of the conservation law of 
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mass, the Bernouhi law for the velocity. So it can be written into the following second-order, 
nonlinear elliptic- hyperbolic equation of mixed type for the velocity potential ip : C M" — )• M^: 

div(p(lVv9p)V(/?) = 0, (2.7) 

where, with a particular choice of the Bernoulli's invariant, the density function p{q'^) is 

piQ') = (l - ^^q') . (2.8) 

We note that equation ()2.7p can be derived formally from the isentropic Euler system. 
Steady isentropic Euler system is obtained from (j2.ip by setting p = Kp^ , where k > is a 
constant, and dropping the third line in (j2.ip . Then equation (j2.7p can be obtained by a formal 
calculation, substituting u = V(p into the isentropic Euler system, and choosing k = ^. 

The second-order nonlinear equation ( 12.71) is 

(2-9) 

strictly elliptic if \S7ip\ > c, and is strictly hyperbolic if \Vip\ < c 



with c = = v^pT-i. The elhptic regions of the equation (j2.7p correspond to the subsonic 

flow, and the hyperbolic regions of (j2.7p to the supersonic flow. 

From the divergence structure of (j2.7p . it is easy to derive the following RH conditions on 
shocks. Suppose that a smooth (n-l)-dimensional surface S divides O into two disjoint subsets 
n^. UipG C°'^(17) is in C'^{n^) n C\n^ U S) then is a a weak solution of ([22]) if and only 
if if satisfies (|2.7p in each of and 



0. (2.10) 



Note that [^p]s = since ip £ C°'^{n). 

Shock S is called transonic if it separates subsonic and supersonic regions. 

Now we discuss radial transonic shocks, i.e. transonic shock solutions of the form 

ip{x) = w{r) where r = |a;|. 

Let < ri < r2 

Q = {ri < \x\ < r2}. 

Then calculation in [CFll pp. 488-489] show that there exists (radial) functions f^{x) = 
w^{\x\ (defined by |CFH eqn (7.4)] with e.g. Rq := ri) such that 

(i) if^ G C~(J7) and satisfy 1^1} in n, 

(ii) is supersonic and 93^ is subsonic in 0. Moreover < max^gf^^ ^.^j --j7-('^) < 

(iii) for any rs € (ri,r2), choosing a constant c = c{rs) so that WQ{rs) = Wq {rs) + c, there 
holds: function 

(fo{x) = m.in{(p'^{x),(pQ (x) + c{rs)) for j; e Q 

is a transonic shock solution of (j2.7p with shock S = {\x\ = r^}, and <po = ip^ in 
0~ := {ri < |x| < r^}; and = tp^ in 17+ := {r^ < |x| < r2}. In particular ipo is 
supersonic in 0" and subsonic in fi+ . 
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Note that dr^p^ > 0. Thus flow enters O through {|x| = ri} and exits through {|x| = r2}. 
By part(iii) above, Vip does not change on the entrance and exit when the position of the 
shock changes within (ri,r2). Since u = V(/J determines the density p by ()2.8|) and pressure 
p = Kp'^ , it follows that parameters of the flow on the entrance and exist of the annulus-shaped 
domain 17 do not determine position r = of the shock. 

Restricting the annulus-shaped domain Q to a cone-shaped (nozzle) domain TV as in Section 
12.41 we see that the the same example shows non-uniqueness of the shock position in the 
cone-shaped nozzle. 

Thus we propose a new potential model, to fix this degeneracy. 

2.3. The non-isentropic potential flow system and transonic shocks. The new model, 
which we call the non-isentropic potential flow system, is obtained by substituting u = V(/J 
into the the full steady compressible Euler system (j2.ip . Unknown functions in the new model 
are the scalar functions {ip,p,p). Note that second line of (j2.ip . representing conservation of 
momentum, consists of n equations. Thus we have overdeterminacy in the case of potential 
model, i.e. we need to retain one equation from the conservation of momentum equations. For 
that, we note that the component in the direction of u of the vector determined by the second 
line of (j2.ip represents conservation of entropy along streamlines (in the smooth flow regions) . 
This motivates us to introduce the following non-isentropic potential flow system: 

div{pSJp>) = 0, (2.11) 

Vv7-dro(pVv7(g)V(/?+p/„) = 0, (2.12) 

div[pv^{]^\vp\^ + = 0- (2.13) 

()2.1ip and ()2.13p represent the conservation of mass and energy respectively, and ()2.12p con- 
cerns the conservation of the entropy along each streamline but it allows the entropy to change 
between streamlines. 

Since (|2.12p is not in a divergence form, we cannot define a shock solution for this new 
model in sense of (12. 3p . Instead, we employ all the R-H conditions (I2.4p - (l2.6p to define a shock 
solution of the new model as follows: 

Deflnition 1 (Shock solutions). {p,ip,p) is a shock solution of the non-isentropic potential 
flow in 0, with a shock S if 

(i) p,pe cO(iF) n c^in^), tp e c^fi^) n c^{n^); 

(ii) (p, '^'P^P)\snTr ~ j^l^'^;P^nn+ ^ ^' 

(iii) {p,ip,p) satisfies (j2.1ip - (j2.13p pointwise in Q^, and (j2.4p - (|2.6p on S with u = \/<p. 
For (/?,(/?, p) e C°°, if p > then (12141 - (|2J3] 1 are equivalent to 

div[{B - ^|V(^|2)^Vv3) = (2.14) 

Vv9-V ^ = (2.15) 

(s-iivv^n^ 

V(/? • V5 = (2.16) 
where the Bernoulli's invariant B is defined by 
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An explicit computation shows that, as a second-order PDE for ()2.14p with (I2.17P is 

strictly elliptic if |V(/?| < c, and is strictly hyperbolic if \Vip\ > c 
with c defined in ()2.2p . Note that this is similar to ()2.9p . 

Definition 2 (A transonic shock solution). A shock solution {p,(f,p), defined in Definition 
is a transonic shock solution if ('i^JdusVlo^cis ^ ^'^"^^n+ns ^ ^ •^^'^ ^^^^ normal Vg on S 
pointing toward , and (ii) |V(^| > c in Q~ , \V(p\ < c in where c is defined in (j2.2p . 

Let Tg be a tangent vector field on S, then, by (j2.5p . there holds 

[p{VifUs){V^-Ts)]s = 0. (2.18) 

Assuming no vacuum state in i.e., p 7^ in 17, ()2.4p and (j2.18p imply either Vy? • i^^ = or 
[V99 • rs]5 = 0. If V(/3 -1/5 = then Definition 13.341 (ii) implies [V(/9 • Ts]s 7^ 0. In that case, 
S becomes a vortex sheet. If [V(/9 • r^Js = holds on S then this implies [<p]s = k for some 
constant k. We are interested in the case of Vi^ ■ Us on S so (j2.18p provides 

[ip]s = (2.19) 

choosing the constant k = without loss of generality. In addition, (j2.5p • Vs provides 

[p{V^-Us)^+p]s = 0. (2.20) 

Remark 2.1. The specific entropy e of ideal polytropic gas satisfies the constitutive relation 

Kexp{e/c^) = 

for constants k and Cy . A straightforward calculation using (j2.4p , ()2.6p , (|2.19p and ()2.20p indi- 
cates that the entropy e increases across a shock S for the non-isentropic potential flow model. 
The entropy jump across a shock plays an essential role for the well-posedness of the problem 
considered in this paper. 

2.4. Radial transonic shock solutions. Fix an open connected set A C S"^^(n > 2) with 
a smooth boundary dA, and define a straight divergent nozzle J\f by 

X 

A/" := {x : ro < \x\ < ri, - — - € A} for < tq < ri < 00. 

\x\ 

Regarding A as a n-1 dimensional submanifold of M"^^, let x' = {x'l,- ■ ■ ,x'^_i) be a coor- 
dinate system in A. Particularly, if A C S"^^, then A can be described by a single smooth 
diffeomorphism(e.g. stereographic projection), say, T i.e., there is A-p C R"~^ with a smooth 
boundary so that A-p can be described by 

A = {T{x') : x' e At}. (2.21) 

Then, for any x G J\f, there exists unique x' satisfying = T[x'). For convenience, we denote 
as 

x=T{r,x') x'StA ifr=|a;|, T(x') = A ^ A. (2.22) 

|x| 

To simplify notations, we denote as x' G A instead of x' ^ hereafter. The coordinates (r, x') 
in (|2.22p is regarded as a spherical coordinate system in A/". 
For later use, we set 

Tent ■■= {x =T {r, x') -.r = ro, x' G A}, Tex ■■= {x =T {r, x') : r = ri,x' e A}, 

^ ex J- 
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To simplify notations, for a function defined for x € A/", if x =t (r, x') then we identify 
^(r, x') with '(/'(x) without any further specification for the rest of paper. Similarly, any function 
/(|^) defined for G A, we identify f{x') with /(y^) if = T{x'). Furthermore, any function 
(fQ varying only by the radial variable we write as '^o{r). Also, in (r, x') coordinates, we 
denote as 

Ar= (ro,ri) X At, M' =Mr^{r <a}, =Mr^{r> a}. (2.23) 



For a fixed constant vector (Pin,Vin,Pin) S 1^1 with Vin > Cin(= a/^t^), let us find a 

V 

transonic shock solution {po,^o,Po) satisfying 

{pO,dr(po,Po) = {pin,Vin,Pin) OU Fe^t, (2.24) 

d^^ipo = onF^ (2.25) 

where is the unit normal on F-y; toward the interior of A/". Since fj2.1ip - (j2.13p are invariant 
under rotations, we expect ipo,(po,Po) to be functions of r only. Then ()2.25p automatically 
holds on F^. For a continuously differentiable radial solution {p,ip,p), ipr 7^ in AA, then 
([2TT]1 - (|2J3]) are equivalent to 

-f(r-W) = 0, 
ar 

d ,1 r.^ dp , ^ 

(j2.26p can be considered as a system of ODEs for {(pr,p) because of 

(7- l)(5o - 2'/'r) 

To find {ipr,p), we rewrite ()2.26p as 



dipr _ 2(re - 1)(7 - l)ipr{Bo - ^(/?^) 
dr ~ (^j + l)r{^^^-Ko) 

dp _ 2(n — l)jip'^p 
d^ ~ ~ (7 + l)r(^2 _ 

with 



(2.28) 
(2.29) 



(2.30) 

Fix E (?^0i ^1)) and let us find a radial transonic shock solution with the shock on {r = r^}. 
For that purpose, we need to solve (|2.28p and (|2.29p in two separate regions of supersonic state 
and subsonic state. 

For the supersonic region M.^^, the initial condition is 

{^r,P){ro) = {Vin,Pin)- (2.31) 

Then (j2.28p . (|2.29p with ()2.3ip have a unique solution (dr^PQ ,Pq) in M^T. For the subsonic 
region A/j;|^, we solve the system of the algebraic equations (j2.4p . (j2.20p and (|2.27p with dr^PQ 
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and Pq on the side of M~ for (dr(p,p){rs)- Then there are two solutions for {dr^p^p){rs) but 
the only admissible solution in sense of Definition [5] is 

dMr,) = (2.32) 
Or^O [rs) 

P{rs) = {poidrfof +Pq - pQKo){rs) =: Psfl{rs)- (2.33) 

So the initial condition for the subsonic region is given by (I2.32j) and (I2.33p . 
Let {drip^,p^) be the solution to <^M), ([2^ with (l2:32]l . dOSjl . 

We claim that {dr(pQ,PQ) is indeed supersonic in 7V"~ while {dr<PQ,PQ) is subsonic in Af^^. 
Plug ^ = into the equation of Bernoulli's law ^{dr(p)^ + (yl^i)p = -^o to obtain 

^^^-c^ = 2±l{^l-Ko). (2.34) 



Denote as c± then, > cf^ and (ITMIl imply ^ > for (^^ = 9.^(^0 thus (5^.(^0 )2 > ci 

for r > ro so {dripQ ,Pq) is supersonic in A/"^". On the other hand, by ()2.32p . dr(pQ{rs) satisfies 
the inequality 

{dr^^f -Cl = I±l— ^(Ko - (dr^of) < (2.35) 

at r = Ts, and this implies < for ipr = drip^ ■ So (dr'PQ)'^ < holds for r > r^. 
Therefore, we get a family of radial transonic shock solutions as follows: 

Definition 3 (Background solutions). For rg G (ro,ri), define 

(Po ,V3o )(^;^^) in[ro,rs] 



{po,po,Po){r;rs) := 



(Po .y'o ,Po )(r;rs) in{rs,ri] 



(7-i){Bo-i(a,.Vo*M)')' 
p~{r;rs):= dr(pQit)dt, Lp^{r;rs) = I dr(pQ {t)dt + ip^ (rs) 



where (9,.(/?o )(r) is i/ie solution to (^TZEtf - l^n^ with (I23D-(I2]MD and ([OT]) respectively. 
We call {po,(po,po){r;rs) the background solution with the transonic shock on {r = rg}. 

Remark 2.2. For a small constant 5 > 0, {dr(p^ ,p^) can be extended to _2s and M^^25 
as solutions to ([2:28]) - ((2:29]) with ([2:32]) - ([233]) and (lOTD respectively where Af^^^g is defined 
by (12:23]) . 

Remark 2.3. For a fixed rg, p^ and p^ monotonically decrease while dr'p^ monotonically 
increases by r. On the other hand, by ([2.35^ . p^ and p^ monotonically increase while dr^PQ 
monotonically decreases. 

Another important property of the background solutions is the monotonicity of the exit 
values of {po.,dr(po,po) depending on the location of shocks r<j. In |YU2j . it is proven that, 
for the full Euler system, p^(ri;rs) and p^(?"i;rs) monotonically decrease with respect to 
while dr'PQ{ri;rs) monotonically increases. If dr'Po{r]rs) ^ in A/", then the non-isentropic 
potential flow model is equivalent to the full Euler system for smooth flow. So the monotonicity 
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properties of the exit data applies to the background solution (poi V'OiPo) of the non-isentropic 
potential flow model so we have: 

Proposition 2.4 ( [YU21 Theorem 2.1]). Let us set pmin =Poi^i]^i) cLnd pmax := Pq^C'^i; ^o)- 
Then Pmin < Pmax holds. Moreover, for any constant Pc € {PmimPmax), there exists unique 
r* G (?^0)^i) so that the background solution {pQ,^o,Po){r;rl) satisfies Pq {ri;r*) = Pc- 

Readers can refer [YU2) for the proof although we would like to point out that the follow- 
ing lemma, which will be used later in this paper, plays an important role for the proof of 
Proposition 12. 4[ 



Lemma 2.5. For any given € {rojri), let (p^, (^J',p^)(r) be the background transonic shock 
solution with the shock on {r = r^} defined in Definition\^ and let us define fiQ by 



(a) Then, there holds /xq > 0. 

(b) Also, we have, for psfl defined in (|2.33p 

rdps,o dp+ (n-l)7Po [(9.9^0-)^ + ^Ko] 

^ dr dr^ ^^^^ - (7 - l)r{Bo - ^dr^^?) ^ ^ ^ 

Proof By l^2B, we have 



d , 



r. 



dr^drifQ 

C{n,^) KoiBo-'^jdrip^f) drift (Bo -'^{drif+f) 

with the constant C{n,j) = "^^""'^^^'^^ > 0. Applying (12.32P to replace 9rvd'(^s) by ^ ^ ^ , 
we obtain 

(2.37) 

dr^dripo (7 + l)rsdriPo (r^) 

By (|2.32|) and (j2.35p . it is easy to see dr{(pQ — ipQ){rs) > 0, and thus (a) holds true for all 

rsG(ro,ri). 

Similarly, by (l^:^ - (ICTD . pTH^ and we obtain ^M^. □ 



-'max ) 1 



2.5. Main theorems. According to Proposition 12.41 for any given constant pc € {pmin,Pr, 
there exists a transonic shock solution whose exit pressure is pc- Our goal is to achieve the 
existence of a transonic shock solution when we smoothly perturb incoming supersonic flow, 
the exit pressure and the nozzle Af. 

Let Af he a nozzle smoothly perturbed from A/" by a diffeomorphism ^ : J\f ^ M", and 
{p, (p,p) be a transonic shock solution of a non-isentropic potential flow in TV with a shock S. 
Then S separates Af into a supersonic region Af~ and a subsonic region Af~^. Let us denote 
{p,(p,p)\j:^± as {p±,'P±,p±)- If the entrance data for {p,(p,p) is given so that the Bernoulli's 

invariant i? is a constant at the entrance of the nozzle Af then, by (j2.4p . (j2.6p and (j2.16p . 



10 MYOUNGJEAN BAE AND MIKHAIL FELDMAN 

(p±i'^±iP±) is a solution to the system of (|2.14|) . (|2.15p and 



(:^^ = ^o (2.38) 



where the constant Bq > is determined by the entrance data. 
On the shock S, (f2lB and (iOSll imply 



[^(V^-P.)^ + ^:^]g = (2.39) 

for a unit normal on S. 

Let us set := ^^ff ( |( V(^- ■ ) ^ + (^Z^ij^, ) , and solve ([23]), UOJ^i and (12:39]) for V(^+-z:;s 
and p+ to obtain 

• = , , (2.40) 

p+= ~P-{V^--Vsf +P-- p-ks on S. (2.41) 

For a transonic shock solution (/), (p,p) in TV, let TV^ := {| V(^| > c} be the supersonic region, 
and TV+ := {|V(^| < c} be the subsonic region where c is defined in ()2.2p by replacing p^p by 
p,^ on the right-hand side. 

Problem 1 (Transonic shock problem in a perturbed nozzle TV = ^'(TV)). Let (p_, ip-^,p-) be 
a supersonic solution upstream with B = Bq at the entrance Tent = ^iXent) for a constant Bq, 
and suppose that V(p- ■ = ^ on = ^'(Fu;) for a unit normal Uw on F^. 

Given an exit pressure function p^x on T^x = ^{T^x), locate a transonic shock S i.e., find a 
function f satisfying 

Kf-=Mr\{\x\ < fi^)}, m = Afn{\x\ > f{^)}, 

\x\ \x\ 
and then find a corresponding subsonic solution (p4-,(^+,p+) downstream so that 

(i) {p+,if+,p+) satisfies (pJil) . (|2T5]) and (12:38]) in TV+ ; ^ 

(ii) (^4- satisfies the slip boundary condition V(^+ ■ = on F^ = dJ\f^ fl F^; 

(iii) i'f+,p+) satisfies (|2l9]) . ([TlO]) and (fTiTT) on S; 

(iv) p+ = pex holds on Tex- 

For convenience, we reformulate Problem [T] for {p,ip,p) = {p,ip,p) o ^ in the unperturbed 
nozzle TV. If {p,'f,p) solves ([2Ti]) . (pT5]) . (p38]) then ip,(p,p) satisfies 

divA{x, L>^', L»93) = (2.42) 

{D'i'^^f{D'^-^)Dip ■ D — = (2.43) 

{Bq - ^\D^-Wip\'^)^ 

hD^-'D^\^+ =Bq (2.44) 

2 (7 - 1)P 

for X € ^'-i(TV±) =: A/'^ with = (dx^if)]^^, = (9x,*j)"j=i, -D*"^ = {D'^yK Here, 
^(x, Z)^', Dip) in (12:42]) is defined by 

A{x,m,r]) = detm {Bq — - \m''^ri\'^)^y^ {m~^)^ m^^r] (2.45) 
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for (x, m,r]) G x B'f'^^\ln) x B^^\o) where we set B^\a) as a closed ball of the radius TZ 
in M.^ with the center at a € M'^. Here, d is a sufficiently small so that ()2.45p is well defined for 
ah m e 

By ()2.19p . the unit normal on S toward Af'^ has the expression of 



with Dy={dy„--- ,dyj. (2.46) 



So the R-H conditions for {p,ip,p) on S := ^^^(5) are 

V9+ = if- (2.47) 

p = p. [{D^-YiD^-^)D^. ■ usf + p_ - p.Ks (2.49) 
with D = {dx^ ) • ■ ■ ) dx„ ) where we set 

Q*i^^^^v) ■■= |J^i(^!'^)| ("^-'f("^-')^ (2-50) 

_ 2(7 -1) ^1 jZ)(v._-(p)|(ZJvI/-i)^j^^-iZJv;_-z., ^, , 7P_ , 

7 + 1 ^2^ |Dxi/-iD(<^__^)[ ^ + (7-l)P-^' ^ ^ ^ 

L>(99_ - (^+) 

= TTv \\- (2.52 

l^(V'- - 95+)| 

We note that, by ()2.47p . z^^ is the unit normal on S toward Af^. Problem [1] is equivalent to: 

Problem 2 (Transonic shock problem in the unperturbed nozzle J\f). Let ^ be a smooth 
diffeomorphism inM. Given a supersonic solution {p-, ip-,p-) = {p-, (^_,p_)o^' inM~ where 
{p-,(p-,p^) is as in ProblemUl and an exit pressure function p^x =Pex°^, locate a transonic 
shock S i.e., find a function f satisfying M~ = {\D'^~^Di.p\ > c} = A/" n {r < f{x')}and 
J\f+ = {\D^~^Dip\ < c} = Af n{r > f{x')} for the (r, x') coordinates defined in (I2.22j) . and 
find a corresponding subsonic solution (p_(_, (^+,p+) in so that 

(i) (p+,v9+,p+) satisfies ([232])- ([231]) in M"^ ; 

(ii) On r+ := n dJ\f~^ , satisfies 

A{x,D^,Dip) ■ = (2.53) 

for a unit normal on T^; 
(in) {if+,p+) satisfies dTTzD - dTlQ] ) on S; 
(iv) p = pex holds on T^x- 

Our claim is that if a background supersonic solution [p^ , (p^ , Pq ) and a constant exit 
pressure £ {Pmin,Pmax) are perturbed sufficiently small, and also if ^' is a small perturbation 
of the identity map, then there exists a corresponding transonic shock solution for a non- 
isentropic potential fiow. 

To study Problem [21 we will use weighted Holder norms. For a bounded connected open set 
17 C M", let r be a closed portion of dQ. For x, y G 0, set 

6x := dist{x, F) 5x,y ■= m\n{5x, Sy). 
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For A; G M, a G (0, 1) and m € Z+, we define 

0<|/3|<m^^^ 

LjC^.r) ._ ;,max(m+a+fc.O) 1-0^^(3^) " D^u{y)\ 

ll„l|(fe>r) ._ II ||(fc,r) I ,(fc,r) 

where we write Z)^ = • • • dx" for a multi-index /3 = (/3i, • • • , /3n) witli /3j € Z+ and |/3| = 
Sj=i /^i- We use tlie notation C^'p^(il) for the set of functions whose || • norm is finite. 

Remark 2.6. Let Q be a connected subset of J\f, and T be a portion of dQ. For T in (|'2.'22p . 
let us set r^^f := {(r, x') : x =t {r,x'),x € Q} and := {(r, x') : x =t {r,x'),x S F}. Then 
there is a constant Cif{N',il,,n,m,a,k) so that, for any u € C™p"(il), there holds 



For a fixed constant data {pin,Vin,Pin) G 1R+ satisfying > ^ and G {Pmin,Pmax), let 
(Po, V^cPo) be the background solution with (po, 9rV'o,Po)(''o) = (/^m, Vin,Pin) and po(?^i) = 
and let Sq = {r = rg} (lAf he the shock of (/Oq, V^OiPo)- Now, we state the main theorems. 

Theorem 1 (Existence). For any given a € (0, 1), there exist constants ai,6,C > depending 
on {pin,Vin,Pin),Pc, ^o, Ti, 7, n, A and a so that whenever < a < a\, if 

(i) a difjeomorphism : AT ^ M" satisfies 

?i := 1^' - Id\2,a,jV < 

(ii) {p^,ip-,p-.) is a solution to ()2.42p - ()2.44p inJ\f~_^_2g, and satisfies the boundary condi- 
tion (j2.53p on F^ fl -N"^ j^25 ^'^'^ ^'^^ estimate 



(iii) a function Pex defined 071 '^ex SatisflGS 



,1 ||(-",9A) ^ 

i/ien i/iere exists a transonic shock solution {p,Lp,p) for a non-isentropic potential flow with a 
shock S satisfying that 

(a) there is a function f : A ^ so that S, J\f~ = {x £ Af : [D'i!^^ Dip\ > c} and 
J\f+= {x eAf : \D^^'^Dip\ < c} are given by 

S = {(r, x')eN':r = f{x'),x' E A}, 

M- = Mn {(r, x'):r < f{x')}, AA+ = n {(r, x') : r > f{x')}, 
and f satisfies 

1/ - rs\i,a,A < C{<^i + + ?3) < Ca; 

(b) {p,ip,p) = {p^,ip-,p-) holds in ; 
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(c) {p^^-,p)\j\f+ satisfies (I2.47p - (I2.49|) . (|2.53|) and the estimate 

||„ I 11,^ _l_ lln ^1+ II ^ r'l'r^ J-r- 

\\P~ P0\h,a,M+ +\\'^~fo\\2,a,M+ + H^' ~ ^0 II l,a,Ar+ < (?1 + ?2 + ?3 j , 

(d) p = pex /loWs on Tex- 

Remark 2.7. /n t/iis paper, we phrase any constant depending on pm, Vin,Pin, Pc, fo, 1^1,1, n, A 
and a as a constant depending on the data unless otherwise specified. 

Theorem 2 (Uniqueness). Under the same assumptions with TheoremUl for any a G (^jl)? 



there exists a constant (72 > depending on the data in sense of Remark \2. 7| so that whenever 
< 0" < (T2, the transonic shock solution in TheoremUlis unique in the class of weak solutions 
satisfying TheoremU\(a)-(d). 

Remark 2.8. By Theorem U^c), Theorem implies that the transonic shock solutions are 
stable under a small and smooth perturbation of the incoming supersonic flow, exit pressure 
and the nozzle. 

2.6. Framework of the proof of Theorem [T] and Theorem [2]. By the decomposition 
(|2.42p - (j2.44p . if we solve ()2.42p for ip, then we can pkig it into the transport equation to solve 
for p, and solve the Bernoulli's law for p. To solve ()2.43p for p, however, the initial condition 
on S can be computed by ()2.48p and ()2.49p . But if we also fix p = p^x on Tex, then the problem 
gets overdetermined. So, we set the framework to prove Theorem [1] and [2] as follows. 



(Step 1) To prove Theorem [H we first solve 

Problem 3 (Free boundary problem for if). Given the incoming supersonic flow {p^,Lp^,p-) 
and a function Vex defined on T^x, find f G C^^^q^j(A) and <p in A/""*" so that (i) ip satisfies 

\D^-^Dip\ < c and ([TIS]) in J\f+ =J\fn{{r,x') : r > f{x')}, (ii) also satisfies (E^TD, dTiS]) . 
([233]) and M{x, Dip) — v^x Tex with the definition of A^(x, Dp) being given later. 

Let us emphasize that Lemma [2. 5 1 is the key ingredient to prove Lemma [3.3l and this lemma 
provides well-posedness of Problem [3l 



(Step 2) In Section [H we solve ()2.43p for p after plugging p in to (|2.43p by the method of 
characteristics in with the initial condition ()2.49p on 5*. However, as we will show later. 
Dp is only in up to the boundary of A/""*" so the standard ODE theory is not sufficient to 
claim solvability of (I2.43p . For that reason, we need a subtle analysis of Dp to overcome this 
difficulty. 



(Step 3) For a fixed (^',99_,p_), define V by ■p(^', (/?_ Wgx) = Pex where pex is the exit 
value of the solution p to (I2.43P where p in Problem[3]is uniquely determined by (^', p-,p-,Vex)- 
In Section m we will show that V{^, p-,p-, •) : Vex ^ Pex is locally invertible in a weak sense 
near the background solution. The weak invertibility is proven by modifying the proof of the 
right inverse mapping theorem for finite dimensional Banach spaces in [SZ] . and this proves 
Theorem [H 



(Step 4) To prove Theorem [21 we need to estimate the difference of two transonic shock 
solutions in a weaker space than the space of the transonic shock solutions in step 3. This 
requires a careful analysis of solutions to elliptic boundary problems with a blow-up of oblique 
boundary conditions at the corners of the nozzle J\f. See Lemma 16.31 
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3. Problem [3l Free Boundary Problem for 99 
To solve Problem O we fix the boundary condition on the exit T^x as 

A{x,In,Dip)Dip ■ 

^ex — Vex 

(3.1) 

where small perturbation of the constant Vc given by 

Vc ■■= A{x,In,D(p^)D(p^ ■ Vex\r,^ = {Bo - ^\drfoiri)\'^)~dr(pQ{ri). 

Here, I'ex is the outward unit normal on T^x so Vex is r = y^. Let ipo,fo,Po) and be as in 
Theorem [TJ 

Proposition 3.1. Under the assumptions of TheoremU^i) and (ii), for any given a G (0, 1), 
there exist constants 0-3, 5,C > depending on the data in sense of Remark \2.T\ so that whenever 
< cj < CT3, if a function v^x defined on T^x satisfies 

?4 := \\Vex - ^c|li,a,A < O", 

then there exists a unique solution ip to (j2.42p . (j2.47p . (|2.48p . (I2.53p . (|3.ip with the shock 
S, and moreover, 

(a) there exists a function / : A — t- so that S, = M Ci {ID"^^^ Dip\ > c} and 
7V+ =N'ri{\D^^'^Dip\ < c} are given by 

S = {(r, x')eM ■.r = f{x'),x' € A}, 

A/-- = A/- n {(r, x'):r< f{x')}, = M D {(r, x') : r > f{x')}, 

and f satisfies 

\f - rs\i,a,A < C{<;i + ?2 + ?4) < Ca, 

(b) (f = (f- holds in J\f~ , 

(c) 93 satisfies the estimate 

y - ^oi'aM^'^^^ ^ ^(^1 + ^2 + < Ca. (3.2) 

We prove Proposition 13.11 by the implicit mapping theorem. For that, we first derive the 
equation and the boundary conditions for ip := — ip^. 

3.1. Nonlinear boundary problem for ip = ip — ip'^ . If is a solution in Proposition 13.11 
then ip := (p — ip'^ satisfies 



dk{ajk\ 


[x, D-ilj)dj4)) = 


dkFk{x 


,D^,Di,) 


in A/'+, 


(3.3) 


{ajk{x, 


D'4))djilj) ■ Vu, 


= [Aix, 


In,D^)- A{x,D-^,Dp)] 


on r+. 


(3.4) 


{ajk{x, 


DTp)djlp) ■ Vex 


— Vex ~ 


- Vc 


on Tex 


(3.5) 



with r+ := n aAA+ and 

ajk{x,ri)= dr,^^Aj{xJn,Dp^ + ar])da (3.6) 
Jo 

Fk{x,m,r]) = -Ak{x,m,Dp'^ + rf) + Ak{x,In,Dip'^ + rf) (3.7) 

where A = {Ai, ■ ■ ■ , An) is defined in (j2.45p . The boundary condition for ip on S needs to be 
computed more carefully. 
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Remark 3.2. In this paper, we fix the Bernoulli's invariant B as a constant Bq. So, for given 
(^',99_,p_), by (|2.44|) . /)_ is given by 

^ IP- 

P- (7-l)(i?o-i|dM/-lZ)(^_|2)- 

So we regard as factors of perturbations of the nozzle, and the incoming supersonic 

flow. 

Lemma 3.3. Let f and (p be as in Proposition \3.1\ and set ip- := (p- — cp^ . Then there are 
two functions fi = fif, gi = gi{x,D'^ ^p^^ip-^Dijj-^Dil)) in sense of Remark \3.S\ and a vector 
valued function hi = bi{x,Dip-,D'tl') so that ip satisfies 

bi ■ Dip — fiftp = gi on S. (3-8) 

Moreover, and gi satisfy, for /iq > defined in Lemma \2. 5t 

II/"/ - /"o||l,a,A < C{c,i +(,2+ (.a), (3.9) 
11^1 - '^^llJXr^ ^ ^(-^l + ^2 + ?4), (3.10) 

hit:r^ <c{,i+<^2). (3.11) 

So if as (in Proposition VJ. 1\) is sufficiently small depending on the data, then there is a constant 
A > satisfying 

fJ-f > — > 0, and bi ■ Us > X on S (3.12) 
where the constant A also only depends on the data as well. 
Proof. Step 1. In general, at {f{x'),x') for x' £ A, tp = ip — p'^ satisfies 

drip = drP - + {^^^ - drP^) (3.13) 
OrPo OrPo 

where Kq is given by (pIHOjl . By (f2:32]l and (f2:T7l) . (fXT3l) is equivalent to 

Dip ■ f - iifip = h (3.14) 

with 

, !iUiri--dr^t)irs + tU{x')-rs))dt 

ljij(^x ) := , (3.15) 

Si driPo - Pt)i'^s+t{f{x') -rs))dt 

h = Dp ■ f ^ fJ-fif- - fo) ='■ h* — fJ'fip- (3.16) 

Dpo -r 

where we set -0- =: P- ~ ■ 

Suppose that = Id, {p^,p^) = {p^ ,p^) and p in Proposition EH is a radial function with 
a shock on S = {r = /g} for a constant /o G (ro,ri). Then, (j2.32|) implies /i = in (j3.14|) . 
Then the boundary condition for ip on S is 

drip - /i/oV' = on S = {r = fo}. (3.17) 

Step 2. (|3.17|) shows how to formulate a boundary condition for p in general. 
Since J^°_ — drP^ is smooth, by Proposition 13. H we have 

II/"/ - /"o||l,a,A < C\\f - rs||l,a,A < C'||V' - V'-lll.a.S < C{c;i + ?2 + q)- 
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()3.9p is verified. By Lemma 12.51 no is a positive constant. 

We note that h* in (|3.16|) depends on D^,p-, Dip- in sense of Remark l3.2l and also depends 
on Dip. So we should carefully decompose h* in the form of /i* = /3 • Dip + g so that a weighted 
Holder norm of g does not depend on ||V'|l2 a '^^ sufficiently small. For h* defined 

in (I3J6]1 . by (H^H]), we set h* = ai - 02 with' 

ai = 1)99 • r - QDip ■ Vg, 02 = ^ , - ^r^^ on 5*. (3.18) 

Di^Pq ■ r QDip ■ Vs 

We note that, by (j2.47p . i/^ — r is written as 

Vs-r = v{Dip_,Dip) - 1/(0,0) = ii{Dip.,Dip)Dip_ + j2{Dip_,Dip)Dip (3.19) 
where we set 

ji(e,r?):=/ D^u{t^,r,)dt , j^i^^n) := [ D^u{0,tv)dt (3.21) 
JO JO 

for ^, € K". Then, we can rewrite ai as 

oi = D(f ■ (r - z/5) + (/ - Q)D(f ■ Us 

= -[g{DiP_,DiP)]D{ip+ +ip).DiP + hi := Pi- Dip + hi 
with 

hi{D^,DiP^,DiP) = -Dip-[ji{DiP^,DiP)DiP^] - {Q - I)D^ ■ u^. 

Here, all the quantities are evaluated on {f{x'),x') € S. 

To rewrite 02 for our purpose, let us consider Kq — Kg first. 

Ko-Ks = l^i iD^,.rf-iQD^_..g)^) + - ^) (3.23) 



(3.22) 



||(^ ^)|5|j[ a A^^ independent of (/? but should be treated carefully. 

= \Dipo\\l - (f • Ugf) + [{D^Q ■ Vsf - iQD^_ ■ u^)^] =: h + h 

It is easy to see that a Holder norm of I2 is uniformly bounded independent of ip. Let us 
consider Ii now. By (|2.52p . a simple calculation provides 

where a Holder norm of Is is independent of ip if is sufficiently small. By ()3.18p . (I3.23P and 
(13211), a2 is 

^ {-f-l)driPo{DlP-drlPf) , [iJ(jJV;_,iJV^)]K,f ^ 

- ^ (7 + l)|I?(^--V^)P + QI?^- • ^ • + (3.25) 

= : P2-DiP + h2 
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with 



7 + 1 (7 + 1)5^(^0 p- 

{Q - I)Dip_-Vs + Djj^Us , i-Tfr,i / ^ 



(3.26) 



QDip^ ■ Us dr<fQ 

By ()3.22p and (I3.26p . we can write h* = ai — a2 as 

h* = (/3i - p2) ■ Dip + {hi - h2) (3.27) 

for /3i,2 = l3i^2iDip-,Dip),hi^2 = /ii,2(^^, ^'V'-, -D^) defined in (|3:22D and i^M)- From this, 
we define bi , (71 in (j3.8p by 

61 =f-/3i +/32, c/i = /ii - /i2 -p/V'-- (3.28) 
By the definition of gi in (]3.28p and Proposition 13. 1^ we have 

\\gi{D^,p^,Di^.,i^.,Dij)\\[-^f''^ <C{l + a)ic^i+<i2) 

so p. lip holds true if we choose 0-3 satisfying < 0-3 < 1. 

It remains to verify (j3.10p . From the definition of u{(,,ri), one can check 

so we have j2(0,0)r = where j2 is defined in (|3:2TD then (13:22]) and (13:25]) imply /3i(0,0) = 
/3i(0,0) = 0. As vector valued functions of ^,77 G M", /3i and P2 are smooth for ^, ?y G -B^(O) 
for a constant (i > sufficiently small. So if 173 is sufficiently small, then we obtain 

||/3i,2(D^-,i?^)|lS7„"r^ < C(ft + ?2 + q), (3.29) 

and this implies 

\MDi;.,Di;) - ^.H^X^ < C(,i + ,2 + ?4) + - rllS;^^^ 
By ()2.52p and Proposition 13.1] Vs — satisfies 

therefore (j3.10p is verified provided that we choose (T3 in Proposition 13.11 sufficiently small 
depending on the data. □ 



3.2. Proof of Proposition 13.11 Fix a G (0, 1). For a constant i? > 0, keeping Remark [3: 
in mind, we define neighborhoods of Id, {ipQ ,Pq), Vc, r^. 

f?(Ji)(/d) := G C2'"(Ar,M") : \^ - Id\2,a,X < R}, 



B^i'\v,) := {ve. G Cl:_l^eA)W : \\ve. - v^Wi^^f^^ < R}, 



(3.30) 
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For a sufficiently small R > 0, B^^\ld,ipQ ,Pq ,Vc) is a set of small perturbations of the 
nozzle, the incoming supersonic flow and the exit boundary condition for ip. For two constants 
C*,a* > 0, define 5 : Bi\\ld,^^ ,v,) x 6gi^,(r,) ^ C^'^^.^ (A) by 

^-,P-,Vex, f) := {^^ - V'o - i'){f{oc'),x') (3.31) 
for x' € A where ^|J satisfies the following conditions: setting 

M+ :=Afn{r> fix')}, Sf := {r = f{x'),x' G A}, 

(i) ip satisfies 

for a constant M > to be determined later; 

(ii) ip solves the nonlinear boundary problem: 

N+{^) := dk{ajk{x, D^)d,ilj) = ^ D*, D^) in 

j,k k 

= gi{x, D'^jP^jip^, Dijj^, Dip) on Sj 
M^iiP) := [ajkix, Di,)]DiP ■ = g2{x, D^, D^j) on T^j := dMf n F^ 



(3.33) 



with 



M^(V') := [ajk{x,Dil))]Dil) ■ v^x = gsivex) on Fg 



52(x,L"I',DV') = -[A{x,D^,Dip+ + Dtp) - A{x,In,Dip+ + Di;)], 

g2.{Vex) = Vex " Vc 



(3.34) 



where ^ in is defined as in For ct* > sufficiently small, if M in ([3:32]) is 

appropriately chosen depending on the data so that ()3.33p is a uniformly elliptic boundary 
problem with oblique boundary conditions, then Z is well defined. More precisely, we have the 
following lemma. 

Lemma 3.4 (Well-posedness ()3.33p ). Fix a E (0,1) and C* > 0. Then there is a constant 
M, o"" depending on the data in sense of Remark \2. 7| and also depending on C* so that the 
fallowings hold: 

(i) for any given {^,ip.,p^,Vex) S b''^^ {Id,ip^ ,Pq ,Vc), f £ 6^i^j(rs), p.33p is a uni- 
formly elliptic boundary problem, and it has a unique solution ip satisfying ()3.32p . 

(ii) i/ o"" is sufficiently small depending on the data and C* , then there is a constant C 
depending only on the data so that the solution of (I3.33|) satisfies 

ii^ii2;xf'"^^^('^i+^2+?4) 

where ?i , ?2 , ?4 are defined in Theorem \2.4\ and Proposition \3.1i 

We will prove Lemma 13.41 bv the contraction mapping principle. For that purpose, we first 
need to study a linear boundary problem related to (|3.33p . Let u* > be a small constant to 
be determined later. For a given constant C* > 0, fix (^, fg^^) E B^^} {Id,ip'^ ,Pq ,Vc), 
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(2) — 

f € Bjj,^:,{rs), and let us set ijj- = ip- — (p^ . Reduce a* sufficiently small depending on the 
data so that D'^ is invertible in J\f, and set 

}Cj{M) := {0 G C^'^Mf ) : M^'^J^f^^ < Ma*}. 
For a fixed 4> G /C/(M), consider the following linear problem in A/^": 

Y,dkiajk{x,D<j))d,u) = ^dkFkix,D^,Dcl)) in A(+, 

bi{D4>,(t>) ■ Du- fifu = gi{x,D'i/,p_,i{j_,Dtl;^,D(j)) on 5/, (3.35) 

{ajk{x, D(f>)dju) ■ Uu, = g2{x, D^, D(f)) on dN'j n T^, 

(ajkix, D(j))dju) ■ Vex = gsivex) on 

with tjj- := 99_ — . The definitions of ajk,bi, fif, gi, §2, gs are given in (j3.6p . (j3.15p . (|3.28p . 
and 

Lemma 3.5 (Well posedness of (|3.35p ). For any given a G (0, 1), and M, C* > is a* > is 

chosen sufficiently small depending only on the data, C* and M , then there exists a constant 
C depending on the data and a so that, for any given f G -6^4* (r,) and G /C/(M), (1335]) 
has a unique solution u^'^'^ G C(^l" _q r )^^t'^' Moreover u'^'^^ satisfies 

\W^'^^i~^;^f''^<C{^i+^2 + <^,). (3.36) 

To prove Lemma [3.51 we basically repeat the proof of |CF3l Theorem 3.1]. A key ingredient 
of the proof in [CF3J is the method of continuity and the divergence structure of a linear 
boundary problem. In our case, however, (|3.35p is not in a divergence form because of the 
derivative boundary condition on Sj. For that reason, we need Lemma 13.81 so that (j3.35p 
can be treated as a divergence problem. Before we proceed further, we first note two obvious 
lemmas for the later use. 

Lemma 3.6. For any C*, M > 0, and for f G B^^}^,{rs), (j) G /C/(M), there hold 

\\ajk{x,Dcl)) - ajk{x,0)\t''f.X^ < CMa*, (3.37) 

\MD<p,<p)-hiO,0)\\[-;;'^J^ < CMa\ (3.38) 
II/"/ - /"o||i,Q,A < C\\f - rs||i,a,A, (3.39) 



||F(x,I)^')|||^^-^^-) + ||5i(x,DvI/,p_,V._,I)V'-,^'^)ll!:X^ 

{-a,9A) 



+ \\g2{x,D^,D<l,)\Ur^^^ + ||<73(^^ex)|lUT ^ (^-^^^ 
where the constant C in (j3.37p - (|3.40p only depends on the data in sense of Remark \2. 7\ Thus, 
for any C*,M > 0, if a* is chosen sufficiently small depending on the data, C* and M, then 
the equation of ()3.35p is uniformly elliptic in Mj' . 

Proof. ()3.37p - (|3.40p are obvious by the definitions. 
For any given ( G M", aij{x,0) satisfies 

aijix,0)QCj > |^^(^o - ^\D^+\')'^{Ko - \Dipt\')\C\'. (3.41) 
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So, by p.34p . ()2.35p and Remark 12.21 for any C* , M > 0, one can choose a* depending on the 
data, C* and M so that (I3.35P is uniformly elhptic in Afj~ . □ 

Lemma 3.7. For any C*,M > 0, let vj be the inward unit normal on Sf, and bi be defined 
in (|3.28p . Then, for any given/ € S^,^*(rs), 4> G K,f{M), we have 

biiD(f>,(f>)-Uf>^ on Sf, (3.42) 

/i/>^^0>0 on Sf for any f e Bgl^rs) (3.43) 

where fiQ is defined in Lemma \2.5\ if a* is chosen sufficiently small depending on the data, C* 
and M. 

Proof. For i/ defined in ()3.20p . we have 
bi{DcP,cP)-Uf 

= ibi{Dc^, <P) - 6i(0, 0)) • uf + 6i(0, 0) • (z// - u{0, 0)) + 6i(0, 0) • u{0, 0) 
> -C{M + C*)a* + 1 

for a constant C depending only on the data. By (I3.38p . for any C*,M > 0, if we choose a* 
to satisfy -C{M + C*)a* + 1 > ^, then (f3:i2D is obtained. Reducing a* if necessary, dOOjl 
provides ([3:i3l) . □ 

By Lemma 13.61 and 13.71 for any C*,M > 0, choosing a* sufficiently small depending on the 
data, C* and M, (j3.3p is a uniformly elliptic linear boundary problem with oblique boundary 
conditions on the boundary. 

Proof of Lemma 13.51 

Step 1. Suppose that u solves (I3.35P and then we first verify (I3.36P by following the method 
of the proof for [CF31 Theorem 3.1]. First, we apply the change of variables x i— t- y = <I>(x) 
with defined by 

\y\=rs+ "''//A lxl-fix')), ^ = n' (3-44) 
n-f{x') \y\ \x\ 

This change of variables transforms Aff' to M^^, and Sf = {r = f{x')} n Mf' to = {r = 
rs} n 7V"^ in the spherical coordinates in (j2.22p . Then w{y) = no solves a linear 

boundary problem of the same structure as (|3.35p . For simplicity, we write the the equation 
and the boundary conditions for w as 

du{ajk{y,f,Dcl))djVo) = dkFk{y,f,D^,D^) in M+, 
bi{y,f,(l),D(l})-Dw-HfW = gi{y,f,D'^,p^,^^,DTp_,D(p) on Sq 
{ajk{yJ,D(/))djw) ■ = g2{D'^,f,D(p) on F^^^^, 
iajkiy,f,D(j))djw) ■ Uex = hivex) on Tex- 

with Tw,ri, = '^wf^dM^^ and ip- = ip^—ip^, where cijk-, bi,Fk,gm are obtained from ajk, bi,Fi^,gm 
in ()3.35p through the change of variables x y = ^(x) for j,k = ,n, and m = 

1,2,3. Next, we rewrite p.45p as a conormal boundary problem with coefficients ajk{y,0) 
and 6i(y,0) = r for the equation and the boundary conditions of w. 
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We remind that, by ()3.6p . the matrix [ajk{x,0)] is 



2-7 
I 7-1 



K-.(y,0)] = {Bo - ^\D^+\')^In - ^-^^^^^^^^^{D^tr D^t- (3-46) 

Lemma 3.8 (Boundary condition for w on Sq). Suppose that w solves (|3.45jl . and let M, C*,a* 
be as in Lemma WTR For any given a G (0, 1), / € B(jl^t.{rs), (f> G K,f{M), there is a constant 
C depending only on the data in sense of Remark \2. 7\so that w satisfies 



{ajk{y,0)djw) ■ Us = hi on Sq (3.47) 
for a function hi G ^^(^1^ gA)(^)' '^'^^ satisfies 

ii/^iiiS:X^ < cm+cna^w\\[-^-^f-^ + \m[-;^;^t^ + ii^iiitX^). (3.48) 

Proof. The unit normail on Sq toward J\f^^ is r. Then, by (j3.46p and the boundary condition 
of w on 5*0 in (j3.45p . there holds 

(a,fc(y,0)a,n) ■ = {Bo - ^|D^+P)^(1 - _ ^^ll^^^l^j^^+^2^ )drW 

= ; —[{r-bi) ■ Dvu + HfW + gi\ =: hi 

{^-l){Bo-^2\M\')~ 

By p.38p and the definition of bi, (I3.48P is easily obtained. □ 
Back to the proof of Lemma [3. 51 using Lemma [321 ^ solves the conomal boundaru problem 
dk{ajk{y,0)djw) = dk[Fk{y, f ^Dcj)) + Sojkdjw] =: dkGk in 7V"+, 
{ajk{y,0)dju) ■ys = hi on So, 

dk{ajk{y,0)djw) • f ^ = 52 + {5djkdjw) ■ Uyj =: /i2 on Tyj^rs, 
dk{ajk{y,0)djvu) ■ Uex = h + {Sdjkdjw) ■ Uex =■ h^ on Tex 
with ddjk := ajk{y,0) — djk{y, f, Dcp). We note that, by (|3.37p . 5djk satisfies 

¥a,k\f-:ij^ < C{m[-^J,f-^ + 11/ - r4;^'r'''^) (3.50) 

for a constant C depending on the data in sense of Remark 12.71 

Step 2. Using a weak formulation of the equation for w in ()3.49p with ()3.46p . one can check 

dk{ajk{y,0)djw) = dkGk 

n-l 

<^ Lq := dr{ki{r)u)i{x')drw) + dx'^{k2{r)bim{x')dx'^w) = dr{G ■ f) +^5^./(G' • x[) 

1=1 

with 

{j + l)r-\Ko-{drip+{r)y) , ^ „_3^^ 1 



(3.51) 



ki{r) := ^ ; ^ ^ "7-2 , Hr) ■= r"-^(i?o - ^{dMY)^ 

2{j-l){Bo-l{dr^Ur)r)— 2 

r'^-V(x') = det 



1 



(3.52) 



d{r, x') 
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for (r, x') is given as in (j2.22p . and every bim(x') is smooth in A. Here we set G = ui[^^^^^]'^G. 

We note that uji is bounded below by a positive constant in A, and the hnear operator Lq 
is uniformly elliptic in ,and bounded by p.4ip and (|2.35p . 

Step 3. To obtain the desired estimate of u in (|3.35p . we use (|3.49p - (j3.52p . Lemma [3^81 
and the interpolation inequality of Holder norms to follow the method of the proof for |CF31 
Theorem 3.1], and prove the following lemma. 

Lemma 3.9. There exists a constant k > depending only on the data in sense of Remark 
17\so that whenever Ma*,C*a* < k, for any f G B^\rs), (j) ^ /C/(M) with a € (0, 1), there 



is a constant C depending on the data so that, if w solves (j3.45p then there holds 

- ^(\'^\o,^^+ + \^\o,aMr+ + l5llo,a,5o + 1^2 lo,a,r„,,), + l^slo.a.FeJ- (3-53) 

The only difference to prove Lemma (|3.9p from |CF31 Theorem3.1] is that, to apply the 
method of reflection as in |CF3t Theorem3.1], we need to extend ajk{y,0) in r-direction by 
reflection on and Tex- For that purpose, we use (I3.52p . Except for that, the proof of Lemma 
21 is basically same as the proof for [CF31 Theorem3.1] so we skip details. 
Once (13.531) is obtained, then the standard scaling provides C^l"_^ p ^ estimate of w. For a 



fixed point yo € A/"rt \ F^^.^^ let 2d := dist{yo,Tw,rs)i and we set a scaled function w^^°\s) near 
2/ = yo by 

(yo)r \ w{yo + ds) - w{yo) 



for s G := {s G Bi{0) : yo + ds £ JV^t \ F^.rJ. By [GTI Theorem 6.29] and Lemma E? 

for each yo ^ A/"rt \ Fu,^.,.^, ^(J^o) satisfies 

- ^^\'^\o,N+ + ?1 + ?2 + ?4) 

where the second inequality above holds true by (|3.1ip . (|3.34p and ()3.40p . 

By (j3.39p and Lemma 12.51 reducing a* if necessary, we can have /i/ > ^fiQ > on A. Then, 
by ()3.54p . Remark 12.61 and [LIl, Corollary 2.5], we conclude that if ii; is a solution to ()3.35p 
then it satisfies (I3.36p . 

Step 4. Finally, we prove the existence of a solution to (I3.45p . Consider the following 
auxiliary problem: 

'dkiajkiy, 0)dju) -u = dkFk in A/'+ 

{ojkiy, 0)dju) ■Vs = gi on Sr, ^2 

{ojkiy, 0)dju) ■Vui = 92 on F^,^^ 

, {ajk{y, 0)dju) ■ Vex = 93 on Fg 



ex 



for F G C"(A/'rt), 511 93 £ C"(A) and g2 G C"(F^_rJ- (|3.55p has a unique weak solution 
because the functional 

I[u\ = \ I ajk{y,0)djudku + u^ -2F ■ Du+ [ {F ■ Uout + g)u 

has a minimizer over VF^'^(7V^) with the outward unit normal Uout of 5A/j^ where we write 

9 ■= 9iXSo + 52Xr„.r, - 93Xr,, 
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with XV defined by 

Xv{v) = 1 for y E P, Xv{y) = otherwise for V C M"". 

For F G C^f^j- ^{M^J, 91,93 E C^f^^Q^^{A), the weak solution u also satisfies (|3.36p . By 
the method of continuity, we conclude that the linear boundary problem (j3.35p has a solution 
u € C'(l"_„p ^{Afj^). By the comparison principle and [LIU Corollary 2.5], the solution to 
()3.35p is unique. □ 

Remark 3.10. According to the definition of the laplacian A^n-i on a sphere we note 

where we replace — Agn-i by A^n-i in the widely used convention (e.g. |BR] ) o/ A^n-i. 
Lemma 13.41 is a direct result from Lemma 13.51 

Proof of Lemma \3.4\ For M,C* > 0, let a* be as in Lemma [331 Fix (^', p_, Wex) E 

Po ^"Po^Po '^c) and / G S^]^,(r^), and define a nonlinear mapping by 

where u^''^^ is a unique solution to (j3.35p satisfying (|3.36p . By Lemma 13.51 if we choose M as 

M = 6C (3.56) 

for the constant C in (j3.36p . then maps /Cj(M) into itself which is a Banach space. 

For A; = 1,2, let us set = Gfi'pk)- By subtracting the boundary problem (j3.45p for W2 
from the boundary problem (j3.45p for wi, one can directly show that there is a constant C 
depending on the data so that there holds 

WQA^i) - Gf{^2)\\[~^J^f-^ < Ca*U, - 0211^;'-;^"^ (3.57) 

So if we reduce a* > satisfying Lemma ESI (I3.42|) . (13.430 and Ca* < 1, then the contraction 
mapping principle applies to Therefore Gf has a unique fixed point, say, ip. Obviously, -0 
is a solution to (j3.33p and it also satisfies Lemma 13.4( 11). From this, we choose CT« = a*. We 
note that the choice of o"" depends on the data( in sense of Remark 12. 7p and C* . □ 

By Lemma [331 3) defined in (j3.3ip . is well-defined when we replace a* by fj" in (j3.3ip . 
For any given (^', (/?_,_p_, Vex) E B^^f\ld,ipQ ,Pq ,Vc) and / E 0^i^„(rs), let ip-^ be the fixed 
point oi Qf. Suppose that /* G 13^i^^{rs) satisfies 

y{^,^.,p.,Ve.J*) = {^--^^-i^^'){f*ix'),x') = in A. 
Then 

^'*:=l''l "-^C (3.58) 

is a transonic shock solution of Problem [2l with a shock S = {r = f*{x'), x'}r\J\f. Moreover, ip-^* 
satisfies Proposition I3.ir a)-(c). Therefore, if there is a mapping & : q = (^', "Uex) ^ f 

so that ^{q,&{q)) = holds, then Proposition 13.11 is proven by Lemma [3^ 
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Lemma 3.11. Fix a G (0, 1) and C* > 0. For 2 defined in (|3.31|) . there holds 

^{Id,ipQ ,Pq ,Vc,rs) = 0. (3.59) 
Also, there exists a constant 5^ > depending on the data and C* so that 

(i) 5 is continuously Frechet differentiable in B^^\ld,ipQ ,Pq ,Vc) x B'^lg (r<j), 

(ii) Df2{Id,ip^ ,Pq ,Vc,rs) : C^:^^_^g^~^{A-) C^^l" _c,,aA) (^) invertible. 



with 



Proof. (j3.59p is obvious. 

Fix qo = {^,Lp_,p_,Vex) e B''J^^^{Id,ipQ,pQ,Vc),fe B'^l^y^irs), q= {"^ , (f_,p_,Vex) G 

a 

q\\i := mk.M + y-haM-^^, + \\P-\kaM-^, + W^e.WtT^ = 1, 

?||(-l-a,9A) _ 
/ll2,a,A — 

To prove the Frechet differentiability of 5, we compute 

-Q-J\e=o2{qo + eq,f + ef) = -Q^\e=oiv~ + e^- -^t- ^e){f{x) + ef{x'),x') (3.61) 

where ■0e satisfies ()3.35p in A/"^ - with replacing /, (j), ,p-,Vex by / + e/, V'e, ^ + e^,P- + 

ep-,Vex + £Vex respectively. Since the case of {p^,(p-,p-,Vex) / (0,0,0,0) can be handled 
similarly, we assume {p-,(p-,p-,Vex) = (0,0,0,0) for simplification. 



To compute ^\e=oips{f{x') +ef{x'),x'), we use ^f+^f, given by 

^f+ef ■■ (^'^') ^ ( ''V^ - - + /(^'),^'), (3.62) 

^ n-f-ef 

so that ^£ := -(/^g o {^j^^j)^^ is defined in for every e with |e| < The corresponding 
boundary problem for ip^ is as in (|3.45|) with replacing /, (j), ^,p-,Vex by f+ef, ip£,'i'+e^,p- + 
ep-,Vex + £Vex respectively. 

ajk {x,rj), F{x, m,r]), bi{x,C,v), ai{x,m,p,^,r]), g2{x,m,r]), g3(x,z), defined in I^M), dSZ]), 
p.28p and ()3.34p . are smooth with respect to their variables. So, subtraction the boundary 
problem for iJjq from the boundary problem for ip^, one can show, by the elliptic estimates, that 
converges in C^l^_^ p )(A/j") to u where u solves the following linear boundary problem 

dk{{ajk{x, + ofj^)dju) = dkiaiDf + aaZ?!-) in AT/, 

{bi{D{ip~ -ip7.),D'il)^)+o^'^^)-Du- iifu = a3f + aiDf + a^D^ on Sf, 

n (3) - - (3.63) 

({ajk{x, Dtp ) + o\^)dju) ■ = dQDf + ajD"^ on r^^,j, 

{{ajkix, Dlp^) + 0^^^)djU) ■ Vex = Vex + agDf + agD^ on Tex- 

where o^jl{l = ,4),aj(i = ,9) (whose dimensions can be understood from ()3.63p 

so we only give symbolical expressions on the righthand sides of (j3.63p ) are (vector valued) 
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functions with the dependence of 

ofl=ofl{x,D^,Di;% 

ai = ai{D^, (92-, P-) - {'^0 ,Po),D{lp_ - ipq), Dil^oJ, Df), 
and, by Lemma [331 each o^^ satisfies 

||og(x,D^,Z)V°)|i;;;_j:j^ < Cmf-^-^f-"^ < CaK (3.64) 

Moreover, it is not hard to check that ojl, Oj are smooth with respect to their arguments 
in a small neighborhood of the background state. By the smoothness of o^-^'^ii the 
smallness of o^^ obtained from (|3.64p reducing o"" if necessary, (j3.63p has a unique solution 

By subtracting the boundary problem for u from the boundary problem for , one can 



also check that there exists a constant C so that there holds 11 -(1^^ — ih^) — u\\^ ^ ^ 
or equivalently 



\\i,--^^-eu\\^-^^f-^ <Ce\ (3.65) 

We emphasize that the constant C in (j3.65p is independent of the choice of q and /, but only 
depends on the data in sense of Remark 12.71 

In addition, ojl , Oj continuously depend on their variables in the corresponding norm for each 

component, also the modulus of the continuity is uniform over B^^^{Id, ip^ ,Pq,Vc) x B^^l^^{rs). 
Back to (I3.6ip . we have shown 

||..o5(«o + e^~,/ + ./)(x') 

= f{x')dri^^ -^0- i^0){fix'),x') + (^_ - U)ifix'),x'). 

By (j3.63p . u linearly depends on q and /. So the mapping C defined by 

d 

^ ■ ^ Q^\e=o^{qo + eq,f + ef) 



is a bounded linear mapping from C2'"(AA,]R") x C'^'^^iM;:^^^) x C2'"(7V;^^5) x C^l" ^^^(A) x 

^I'-i-a dA.)^^^ to C^^l" _Q, gy^) (A) where we understand C'^^l".^ ^^^(A) in sense of Remark 12.61 

It is easy to see from (j3.66p . (j3.65p that, for (g,/) € b'^^^ {Id, ,<fQ ,Pq ,Vc) x B^^l^^{rs), 
DZ{q, f) is given by 

D^iq, f) : {q, f) ^ ^\s=o2iqo + eq, f + ef), (3.67) 

and is continuous in {q,f)- This verifies Lemma l3.11l fi). 

Next, we prove Lemma r3.11l fii). By (j3.7p . (j3.28p . (|3.34p . we have, for any tp, 

Fk{x, Id, Dij) = g,{x. Id, 0, DV) = 92{x, Id, D^j) = g:i{v,) = (3.68) 

for k = \, - ■ ■ ,n. Lemma 13.41 (j3.6ip and (j3.68p provide 

Q _ ~ d _ ~ ~ _ 

Q^\e=o^{Id,iPo ,Pq ,Vc,rs + ef) = —\e=o{^Po - ^o)(.^s + ef) = fdri^Pg - iPo){rs)- 
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Thus, we obtain 

Df^{Id,ipQ ,Pf^,Vc,rs) : f ^ fdr{(Po - V^o (3.69) 
By (|2.32|) and (j2.35|) . we have dr{^PQ —^o)irs) > 0. Thus we conclude that Dj^ is invertible 
at {Id,ipQ ,Pq ,Vc,rs). Here, we choose S = ^ where cr" > depends on the data and C*. □ 

Proof of Proposition \3.1[ Let us fix C* = 10 then the choice of a* is fixed accordingly, and 
moreover the choice of cr* only depends on the data in sense Remarkn 12.71 By Lemma 13.111 
and the implicit mapping theorem, there exists a constant < < 5^, depending on the data 

in sense of Remark 12.71 so that there is a unique mapping S : b'^^ {Id, ,Pq ,Vc) 13'^}^*{rs) 
satisfying 

^{q,S{q))=0 for ah q e B^^^ild,^^ ,Pq ,Vc), (3.70) 

and S is continuously Frechet differentiable in (Id, (pQ ,Pq ,Vc). 

By (|2.19|) . Lemma 13.41 and the observation made after the proof of Lemma 13.41 (|3.70p 
implies that = S{q) and c/j-^^'' in (j3.58p satisfy Proposition 13. II (a)-(c). The uniqueness of a 
solution of Problem [3] also follows from the implicit mapping theorem. We note that the choice 
of (73 > only depends on the data because depends on the data and 6^ which depends on 
the data and C* but we fixed C* as C* = 10. □ 

4. Transport Equation for Pressure p 

This section is devoted to step 2 in section \TM 

4.1. The method of characteristics. By the choice of as in Proposition l3.H if is a solution 
to Problem [3l then dr^p is bounded below by a positive constant in M^. So r can be regarded 
as a time-like variable and then we apply the method of characteristics to solve ()2.43p for 
p. According to (|3.2p . however. Dip is not Lipschitz continuous up to dAf'^. So the unique 
solvability of (I2.43P needs to be established first, and we will prove the following proposition. 

Proposition 4.1. Let be as in Proposition \3. 1[ For any given a € (0, 1) there exist constant 
C > and a' £ (0,(73] depending on the data in sense of Remark \2. 7| so that «/?i + ^2 + ?4 < c'; 
and p is a solution of Problem with a shock S = {r = f(x')}, then (j2.43p with the initial 
condition (j2.49p on S has a unique solution p satisfying 

IIP-PollSli^+^<C(ft + ?2 + Q). (4.1) 

We denote as for the subsonic domain {{D'if'^Dpl < c} = {r > f{x')}. We first note 
that p is in up to the wall P^ of J\f~^ away from the corners. 

Lemma 4.2. Given (\I', fex) G ^3^^ (Id, Pq ,Pq ,Vc) satisfying Proposition \3.1\ (i), (ii) 

reducing if necessary, p, the corresponding solution of Problem\^ with a shock S satisfies 

y - ^0 ^^7-^'"^^ < Ci.i +,2 + q) (4.2) 



for ft , ?2 , ?4 in Theorem [7] and Proposition \3.1[ 

To prove Lemma |4.H it suffices to show that the nonlinear boundary problem (|3.33p in the 
fixed domain J\f^ has a unique solution tp so that p = Pq + ip satisfies (|4.2p , and this can 
be proven by a similar method to the proof of Lemma 13.41 except that the scaling argument 
in step 4 needs to be performed for any point yo G AA+ \ (5 U T^x)- So we skip the proof of 
Lemma 14.11 
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Considering r as a time-like variable, it is more convenient to use a changes of variables 
to transform to keeping (|4.2p in the transformed domain . This requires a new 
change of variables: 

G^{r,x') = {{k{ip^ - Lp){r,x') +rs){l - x{r)) + rx{r),x') 

(4.3) 

= : (f (r, j;'), x') =: {f,x') 



if r < r o + 



with a smooth function x{^) satisfying x(r) = < ^ r^-r ^^'^ x'{^) ^ for r > 0. 

I 1 if r > ri 

By ()4.2p , if <;^i + <f2 + Q = o" is small depending on the data in sense of Remark 12.7^ and we 
choose a constant k satisfying 

Co- < min{(v3o - f^){ri), dr{(Po - (Po)irs)}, 

(4.4) 

^ _ n - rs ^ ^ 

^iv^o - )(n)' 

then we have 

a^>)>!l_l£>0 in A/-+_2„ (4.5) 

thus is invertible. 

We choose x' in (j2.22p so that (f, x') forms an orthogonal coordinate system. To specify this 
choice, we write {r,x') as (r, ??) with i!) = (??i, • • • ,'!9„_i) for the rest of paper. And then, we 
apply the change of variables (j4.3p to rewrite (j2.43p in Af~^ as 

n-l 

V ■DE={V ■ f)dfE + ^(F • ^j)d^^E = for (f,??) G (r,,ri] x A (4.6) 

with 

F = J^(D^-i)^(D^-i)L>99, SoG^ = ^ — (4.7) 

(So - il^^-^^'/'P)^ 

where J,^ is the Jacobian matrix associated with the change of variables in (j4.3p . Here, f, -i^j 
are unit vectors pointing the positive directions of r and "Qj respectively. By Proposition 13.11 
Lemma 14.21 and (|4.3p , we have 

Lemma 4.3. Set Vq := J^+DipQ . Then, V defined in ()4.7p satisfies 

/or ?i , ?2 , in Theorem [7] and Proposition \3.1[ Thus, reducing in Proposition I3.il if neces- 
sary, there is a constant uq > depending on the data in sense of Remark \2. 7| so that there 
holds V • f > ujQ > Q in J\f^^ ■ 

The proof is obtained by a straightforward computation so we skip here. 
Lemma 4.4. Let us set := ^ then E^ is a constant. 

Proof. For Vq defined in Lemma HT3l by (j4.8p . Eq satisfies 

Vb • DE+ = {Vo ■ r)drE+ = 0. 
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Since Vb • r > for all r G [rs,ri], and EQ{rs) is a constant by (I2.32p . ()2.33p . therefore we 
conclude E^{r) = E^{rs) = constant for all r G [rs,ri]. □ 

By Lemma 14.31 we can divide the equation (j4.6p by F • f so that the radial speed of all 
characteristics associated with V is 1. From now on, we consider the equation 

"■'^ V 

dfE + TT^^^. E = in N:;y. (4.8) 

We denote 

w'^^^.^,^^^,^^y (4.9) 

V ■ r V ■ r 

The initial condition for E is given from ([239]). Setting Q2 := (-D*"^)^-D^'"^ Eir^-d) must 
satisfy E{rs,-d) = Einti'd) for all ■!? G A where we set 

Eintiif) — — . !M)=(/W,^) • (4.10) 

{Bo - i|i:>4'-i(^')L»99|2)7-i 

Above, /{"&) is the location of the transonic shock for the potential 99 in Proposition 13.11 To 
simplify notation, we write f as r hereafter. By ()4.8p and ()4.10p . Proposition 14.11 is a direct 
result of the following lemma. 

Lemma 4.5. Fix a G (0, 1), and suppose that a vector valued function W = (1, W2, • • • , Wn) 
satisfies W ■ = 

OH '^w^Vs — ^•^r's ^ '^'W fo"^ ^ unit TlOTTJlCll OTl ClTld 
||^||{-a,r^) ||^||{-a,5^ur;^) ^ ^ 

for some constant k > 0. Then the initial value problem 

n-1 

drF + J2 Wj+id^^F = in N+, F{rs,^) = Fo{^) for ^ e A (4.12) 
i=i 

has a unique solution F satisfying 

iii^iiS:;^:^<ciiFoiiS:;r (4.13) 

for a constant C > depending on n,A,rs,ri,a and k. 
Proof of Proposition \4-l\ 

By and Lemma W* = {I, • • • , ^'^"r' ) in gSD satisfies ah the assumptions of 

Lemma 14.51 So, ()4.8p with ()4.10p has a unique solution E. Then the pressure function p is 
given by 

2' 

It remains to verify (|4.ip . By Lemma 14.41 if ^ satisfies (j4.8p then so does E — Eq . Then, by 
(j4.13p . we obtain 

where Eint is given by ()4.10p . 



p=iBo- -\D^-^Dip\^)— ■ {E o G^). 
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Moreover, by ()4.10p . one can explicitly show 

\\Eint - i^o""!! t^r^ < + ?2 + ?4). (4.15) 

This with (I3.2p gives ()4.ip . We choose a' = as by reducing as if necessary to satisfy (14. 2p and 
(lOl) . □ 

Now it remains to prove Lemma |4.5[ Consider the following ODE system 

X(t) = -W(Xit)) for r <t<2r-rs 

^ ' ^ ^ " ~ (4.16) 

X(r) = (r,7?) ^ ^ 

for (r, i9) G U Fga;. The slip boundary condition VF • i/.^, = on r+ implies that any 
characteristic curve X does not penetrate out of N^^. If (j4.16p has a unique solution X{t] r, -i?) 
for any given (r, t?) G 7V+ U Fga; , then we can define a mapping J' : AA^^ U Tex ^ A by 

J : (r,!)) ^ X{2r -rs;r,'d). (4.17) 

If J is differentiable then = i^o o ^ is a solution to ()4.12p . 

4.2. The system of ODEs ()4.16p . ()4.16p has a solution for any given (r, -i?) € (rs,ri] x A by 
the Cauchy-Peano theorem \AM.\ 7.3]. But the standard ODE theory is not sufficient to claim 
the uniqueness of a solution for (j4.16p because W is only in C"(A/rt). To prove the uniqueness, 
we will use the fact that DW is integrable along each solution X to (I4.16p . Denote as X{t; r, ^) 
for a solution to (j4.16p . Hereafter, we call each X(t;r,-d) a characteristic associated with W 
initiated from (r, "i?). 

Remark 4.6. Writing X = {Xi, ■ ■ ■ ,X„), it is easy to check 

Xi{t;r,'d) =2r -t (4.18) 

for any (r, ??) G (r^, ri] x A. So, we restrict (j4.16p for t £ [2r — ri, 2r — r^] because of Xi{2r — 
ri ; r, ■!?) = ri , ATi (2r — ; r, ??) = . 

Lemma 4.7. For any given (r, i?) G (rs,ri] x A, (j4.16p has a unique solution. Moreover there 
is a constant C > depending on A,n and k satisfying 

^dist{{r, 1?), r^) < dist{X{t; r, ^9), T^) < Cdist{{r, ^9), T^) (4.19) 
for t G [2r — ri, 2r — r^]. 



Proof Since is in C"-{Mrs)-, (|4.16p has a solution X{t;r,'d) for any (r, -i?) G (rs,ri] x A. Fix 
r G {rs,ri\. Let X^^) and X^"^^ be solutions to (|4.16p with 

X«(r) = (r,t9i), X^^\r) = {r,^2) (4.20) 
for ■di,'&2 G A. By (|4.16p . we have 

- = 2(iy(x(2)) - Ty(X«)) • (X« - x(2)) 

< C|W(2r - t, OIl-CA)!^^'^ - ^^'f ■ 
In (r, ??)— coordinates, it is easy to see 

(iist((f, 1?), 5*0 U Fex) = minjri — f, f — r^} for {f,'d)^J\f^. 
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So there is a constant C depending on n,Af^^ so that DW satisfies 

\DW{2r -t,'&)\ < Ck{l + (2r - - + (n - 2r + t)-i+") =: m{t) (4.21) 

for any € A and t G [r, 2r — r^] where /c is as in ()4.1ip . Then we apply the method of the 
integrating factor to obtain 

^(e"/XO'it'|x(i)-X(2)[2(t)) <o for tG [r,2r-r,]. (4.22) 

If = ^2 then ([TO]) and ()i:22]) imply X^'^\t) = X^^\t) for all t G [r, 2r - r^]. We have 
accomplished the unique solvability of (j4.16p . 

Since Vl^-z/^ = on F^, any characteristic, initiated from a point on the wall of TVj^, remains 
on the wah i.e., if t?2 G dA, then X^^) lies on T^^. Then (|4.22p implies 

(iist(xW(t),r^) = inf |xW(t) -X(2)(t)| 
i?2e9A 

< C inf |(r,^9i) - (r,t?2)| = Cdist((r, i?i), T^). 

This verifies the second inequality of (14.19P . One can argue similarly to verify the first inequality 
of (ITOjl . □ 

4.3. Regularity of J : (r,i?) ^ X{2r - rs;r,T^) in KTl\i . 

Lemma 4.8. Fix a G (0, 1), and let k be as in Lemma \4-5\ The mapping J defined in (I4.17P 
is continuously differentiable, and it satisfies 

ll^lllJxi < c 

for a constant C > depending on n,A,rs,ri,a and k. 
Proof. First of all, let us assume that J' satisfies 

l^li,o,< ^ ^' (4.23) 
and we estimate [DJ']^^'^Jl^. If J' is differentiable then (|4.16p implies 

(a) drJ{r,'d) = Z(2r - r,; r, i?) - 21^(X(2r - r^; r, i9)), 

(b) a^J = Z(2r-r,;r,i?) 
where Z{t;r, ■!?) is the solution to the ODE system 

Z(t;r,i?) = -DW{X{t;r,^))Z{t;r,i^), 

' (1, 0, • • • , 0) + W{r, 1?) for (a) (4-25) 
{0,6ij,--- ,6n-i,j) for (b). 



(4.24) 



Z{r; r, i?) 



Fix (r,i?), (r',^?') G (r^n] x A with (r,t?) / {r',^'), and consider J := ■ For 

simplicity, we assume rg < r < r' < ri and '& = '&'. Under this assumption, we first estimate 
^'^^^^'^^Z^'l^a^ '^^ because the case of 7^ can be handled similarly. Also, the estimate of 

^'^\{r''d)-(r''t')\" similar and even simpler then the estimate of ^'^'''^^'^^Z^\a^ '^^^ • 
By (|iTT]) and ([03|) . it is easy to check 

2\W(X(2r -r.:r.^)'\ -WiX(2r' -r.:r''d'\)\ 



TRANSONIC SHOCKS IN DIVERGENT NOZZLES 31 

for a constant C depending on n, A, rs,ri,a and k. For the rest of the proof, a constant C 
depends on n, A,rs,ri,a and k unless otherwise specified. 

To estimate — Zi£Lj£)1^ -^g spUt the numerator as 

\Z{2r - r,; r, - Z{2r' - r^; r' , 
< \Z{2r -rs;r,'d) - Z(2r - r^; r', + |Z(2r - r^; r', i?) - Z(2r' - r^; /, ??)|, 

and consider them separately. 

By (j4.2ip and (j4.25p . one can easily check 

\Z{2r -rs;r',T}) - Z(2r' - r,; r', < C(r' -r)". (4.27) 

Let us denote as Zi(t) and Z2{t) for Z{t;r,'d) and Z{t;r',i^) respectively. Then, by (j4.25p . 
z := \Zi — Z^^ satisfies 

z < C{\DW{X{t;r,{^))\z + \DW{X{t-r,^)) - DW{X{t;r',{^))\z^/^) (4.28) 

on the interval Ij-y ■= [2r' — ri,2r' — r^] n [2r — ri,2r — r^]. By (j4.23p . we have 

\DW{X{t;r,{^)) - DW {X {t; r' , ^))\ 

< [ _ , \DWiX{t;r,^)) - DWiX{t;r',m 
- ' ' \X{t;r,^) -X{t;r',{^)\<^ 

Set d := min{(izst((r, t?), r^,), dist((r', t?), r^)} and 

, . \DW{X{t;r,i3))-DW{X{t;r',^))\ 
\X{t;r,^) - X{t;r',^)\° 

For a constant P G (0,1), we estimate and 2:1-'^ using two different weighted Holder 
norms of W. By ()4.19p . for any t € Ir,r', we have 

< Cd^'^(||VF||S""J:V < Ckl^d-^. (4.29) 

Next, set (i*(t) := min{dist(X(t; r, i?), 5o U Fe^:), dist(X(t; r', t?), 5o U Fe^)}. Then (i*(t) can be 
expressed as 

d^{t) = min{2r — t — r^, ri — 2r' + t} for t € /r,r'- 
So (d*(t))~^"^^ is bounded by m^:{t) which is integrable over Iry, then we have 

Z^-^it) < Cm,(t)(||PF||5-"_;^^^=^V-^ < Ck^-^m,{t). (4.30) 

Back to (jMB), by (g^l]), (05I1 . K29h . ()i30]l and the method of integrating factor, we obtain 

l(e-/Xt')'^t'^) < (/ _ rrckd-^m,z'/' 
at 

=^ z{t) < C((r' - rf + (r - r)°d-^ sup z^/^) for all t G 1^,^/ 
^ sup z <C{{r' -rf + d-^'^{r' -rf") 



which provides 

\Z{2r-r,-r,{})-Z{2r-Ts-y,n < (43^) 
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for any /3 G (0,1) and a constant C depending on n,Af^^,k and /3. The case of "i? 7^ can 
be handled in the same fashion. The estimate of ^'^^'J^r'I)-{r''^')\"'^ ^ even simpler then the 

estimate of ^'^'''J^r'^)-^'''^^) ] '^^ ®° ^^iP details and we conclude that, for any constant 
/3 € (0, 1), there is a constant C so that there holds 



By means of similar arguments, one can directly prove G C^(A/^). □ 
Remark 4.9. J is in C^^^°_^^^p )(-^rt) f'^^ ^''^V P ^ (0' 
Now, we are ready to prove Lemma 14.51 



Proof of Lemma \4T5\ By Lemma SSI F := Fq o J is a solution to (jil^ . By (fiT9]) . we 
easily get the estimate ||-F||| "^^+^ < C'||-Fo|li cT^^^ so it is remains to estimate Fix 
(r, 7?) and (r',??') in (rs,ri] x A, and assume (r, ??) 7^ (r',??'). Let us set 

d := ^ min{(i«st((r, {^),T^),dist{{r' , ^'),T^)} > 0. 

Setting [C]i ■■= Cib) - C{a), we have 

\DF{r,'&) -DF{r',i9')\ 

<\[Djp},^\\DFo{J{r,m + \DJ{r',mD^^ 
By (|4.19|) and Lemma 14.81 we obtain 

\DJ{r',§')\\\DF,oJ]^f)^,^\ <Cd''Kj Kp,{Kj\{r,^) - (r',^')l)" 

with Kj = \J\^^j^+ and = ll^o|lS7„°f "^^^ Thus F = Fq o J satisfies (|iT3]) . 

The uniqueness of a solution can be checked easily. If F^^') and F^^^ are solutions to (j4.12p 
then F^^^ — F^^-* is identically along every characteristic X(t; r, -d) associated with W, and so 
= holds in A/'+. □ 

' s 

5. Proof of Theorem [1} Existence 

To prove Theorem [U we prove a weak implicit mapping theorem for infinite dimensional 
Banach spaces. 

5.1. Weak implicit mapping theorem. Fix a constant R > and set 

B^^\ho) := {h € cIlIsa)W --Wh- hoWtT^ ^ ^i' (^-l) 
For B^J\ld,ip^,pQ,Vc) in (f330D . define P : sl^^d )^^?g)(p,)by 

V :{^,^-,P-,Vea:)^p\T,, (5.2) 

where p is the solution to (j2.43p with ip in Proposition 13.11 for the given {^,(p-,p-,Vex)- By 
Proposition 13.11 and 14.11 there exists positive constants and C so that whenever < a < (T3, 
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"P, defined in ()5.2p . is well defined. To prove theorem [H we need to show V ^{pex) 7^ for any 
Pex G (Pc) for a sufficiently small constant a > 0. So, we need the following lemma. 

Proposition 5.1 (Weak Implicit Mapping Theorem). Let Ci and C2 be Banach spaces com- 
pactly imbedded in Banach spaces !Bi and ^2 respectively. Also, suppose that there are two 
Banach spaces £3,^33 with £3 C QS3. For a point {xo,yo) € x C2, suppose that a mapping 
T satisfies the followings: 

(i) J- maps a small neighborhood of (xo,yo) ^1 ^ ^2 to ^3, and maps a neighborhood 
of {xo,yo) in ^1 X to €3 with T{xo,yo) = 0, 

(ii) whenever a sequence {{xk,yk)} C Ci x C2 near (xq, yo) converges to (x,,, y*) in 5Si x *B2, 
the sequence {J^{xk,yk)} converges to F{x*,y^) in ^83, 

(iii) J-, as a mapping from 5Si x ^2 to 5S3 also as a mapping from (Ti x (^2 to C3, is Frechet 
differentiable at {xQ,yQ), 

(iv) the partial Frechet derivative Dj.F{xQ,yQ), as a mapping from 5Si x 5S2 to 5S3 also as 
a mapping from x ^2 to 6^3, zs invertible. 

Then there is a small neighborhood ^^2(^0) of yo in ^2 so that, for any given y G ^^2(^0); ^/lere 
exisis X* = satisfying T{x*{y),y) = 0. 

In |SZj . the author used the Brouwer fixed point theorem to prove the right inverse function 
theorem for finite dimensional Banach spaces without assuming continuous differentiability of 
a mapping. As an analogy, we apply the Schauder fixed point theorem to prove Proposition 
15.11 The detailed proof is given in Appendix lAl 

To apply Proposition 15. II to V, we need to verify that V is continuous, Frechet differentiable 
at {Id,ipQ ,Pq ,Vc) and a partial Frechet derivative of V at the point is invertible. 

5.2. Frechet differentiability of 7^ at (^q = (Id,ipQ ,Pq ,Vc)- To prove Frechet differentiabil- 
ity of V at (q, it sufhces to consider how to compute the partial Frechet derivative of V with 
respect to Vex at Co because the other partial derivatives of V can be obtained similarly. 
For V defined in (15. 2p . let us define TZ and Q by 



n:{^,ip-,p^,Vex) ^ {Bo--\D^-'D^\ 



7-1 



Q : {^,^^,p-,Vex) ^ r- 

n{'i,ip-,P-,Vex) 

If IZ and Q are Frechet differentiable at Co then the Frechet derivative DV of V at C,q is given 
by 

DV^, = n{Co)DQ^, + Q(Co)^7^Co. (5.4) 
To simplify notations, we denote as D^V, D^Q and D^TZ respectively for the partial Frechet 
derivatives of V, Q and TZ at Co with respect to Vex hereafter. 

Lemma 5.2. Fix any a G (0,1), and set 21 := C^^^q^-^{A). As mappings from B^\Co) 
(2) 

to B(j^{pc), Q,TZ and V are Frechet differentiable at Co- Iri particular, the partial Frechet 
derivatives of Q, TZ and "P at Co with respect to are given by 

D^TZ -.we^^aiwe^ (5.5) 

i^^Qw e 21^ a2V'o(rs,-) G 21 (5.6) 

D^Vw G 21 ^ Q(Co)«iw^ + TZ{Co)a2Mrs, •) e 21 (5.7) 
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for any w G 21 with 



ai 



(7 + lK-i(i^o-|5.^o'"P) ^5 
"2 = ; ] : — — o / 4T + " 



(5o-^|Vvp+(r,)|2)— ^dr{iPo -ip^) (7-l)(So-i|5.(/?o P) 
with fiQ defined in Lemma \2.5l and 

Ps,o = Po I^V'o • +Po - Po -^0 (5.9) 

where Kq is defined in (j2.30p . Here, ipQ is a unique solution to an elliptic boundary problem 
specified in the proof. 

Proof. In this proof, the constants Eq and C appeared in various estimates depend on the data 
in sense of Remark 12.71 

Since ()5.7p follows from (jS.Sp and (|5.6p . it suffices to prove (jS.Sp and (|5.6p . 

If a bounded linear mapping L : 21 ^ 21 satisfies 

\\Q{Id,ipQ,p^,Vc + ew) - Q{Id,ipQ,p^,Vc) - ^Lw\\[~"f^^ < o{s) (5.10) 

for any w G 21 with \\w\\[^ a A^^ ~ ^ ^(^) independent of w, then we have L = D^Q. 

To find such L, we fix w G 2t with ll^^lli cTyf^^ ~ compute the Gateaux derivative 

^\e=oQ{Id, ip^,pQ,Vc + ew). 

Fix a sufficiently small constant and, for e G [— EQ) £o] \ {0}, let be a solution indicated 
in Proposition 13.11 with = Id, {ip_,p_) = {lPq ,Pq) and Vex = Vc + £w on Tex in (|3.1P - 

Let us denote the subsonic domain of (pe as M^, the transonic shock as = {r = fd'd)}, 
as and as Gg where <I>/^ and are defined in (|3.44p and (|4.3p . For W^* , defined 
as m let G C(i° r„)(-AArt) be a solution to 

W* • = in AA+ (5.11) 

i?,(r„T?) = ^ -ifei^),^) (5.12) 

(i?0 - ^|V(^,|2) — 

with = Pq{DlPq ■ i/s,e)2 + Pg - p^Kg^e where D = (5,., , ■ ■ ■ ,5ij„_i) and i^s,e is the 
unit normal vector field on toward Af^ , and Kg^ is defined by (j2.5ip with {p^,ip^,p^) = 
(Po j'/'o iPo ) ^nd ^' = Id. By Proposition 14. H the solution Es to (15. lip with (j5.12p is unique, 
and (|5.3p implies Q{Id, ,Pq ,Vc + ew) = Ei;\r^^. Similarly, Q(Co) is the exit value of the 
solution Eq to 

. DE+ = in E^{r,, ^) - P'''^''^ 



{Bo-l\dr^Urs)\'') — 

where Wg* is defined by (j4.9p with (^f, p_, (^_,p_) = (/d, yjfj' , pg , , pg ) , and ps,o is defined 
in (|5.9p . Since is a radial function, we have Wg* = (1, 0, • • • , 0) = r, so E^ is a constant in 
A/j^, and this implies 

^|e=oQ(M</'o ,Po '^c + eu^)= -^\s=oEe{ri,-)= -^\e=o{Ee- Ef^){ri, ■) . (5.13) 
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Note that, by Proposition 13.11 and Lemma 14.21 if eo > is sufficiently small, then whenever 
< |e| < Eq, there exists a unique V'e € C'^'"(7V+) to satisfy 

and moreover, there is a constant C > so that there holds 

' (5.14) 

\\w* - Wo*||^^^;"/™^ + \\w* - < C\e\. 

By (j5.13p . we formally differentiate We ■ D{Ee — E^) = with respect to e, and apply 
(I31]),(jl2|) and (I5H]1 to obtain 

Wo ■ -^\e=o{Ee - E+) = in M+. (5.15) 



Then (|5.15p implies 



— ^e=oQiId,(pQ ,Pq ,Vc + ew) = -^\e=o{Ee - E^){rs) 
d, ( Ps,s pI ^ (^-16) 



(5o-i|D^,|2)7-i (5o-|l^¥'^P)--^^ 



where p^^e is given in (j5.12|) . 

Setpsle,Me),e) ■.= psMeiO),e). By dSHI), dial]) and Ps{e,-,-) is partially differ- 

entiable with respect to e at e = 0. Then this provides 

^\e=oiPs,e-pWeiO),0)=§^\e=OPs{e,rs,e) + {psfl-p^yir,)^\e=ofeiO). (5.17) 

First of all. An explicit calculation using (j2.52p and (j5.2ip implies 

-^\e=0Psie,rs,e) = 0. (5.18) 

Set ipe '■= ipe ° then (j3.44p implies Ve(/£('*?)) ^) = '4' ei^^s ■,'&)■ We note that eV'e solves 
(|3.45p with (j) = +£'(/'£(= y^e) and = ew. Since all the coefficients, including F,gi and g2 
in (|3.45p smoothly depend on cf) and Dcf), ip^ converges to ipo in the norm 
where i/jq is a unique solution to 

9j(ajfc(x,0)afcVo) = in Ar+, Dipo ■ f - ^oV'o = on (5.19) 

{ajk{x,^)dkipo) ■i'w = ^ on Ty,^rs, {ajk{x,Q)dkipo) ■ r = w on Tex, (5.20) 

where ajkjfJ-o are defined in (|3.6p and Lemma [23] Also, by a method similar to the proof of 
Lemma 13.111 we have 

Us o <^-' - M^'^^^f""^ < C\e\. (5.21) 

By ^M, for any e G [-ec^o] \ {0}, we have ((^q " Vo){fem = #e(/e(^?), ^9)- Combining 
this with ([CTTp . i;^:^ and (f^T^ . we obtain 

;^|.=o(p.,e-Po^)(/.(^),^) = i^^'°_"^°7,?\ ^o(r.,g). (5.22) 
Similarly, by using (j5.19p . (j5.2ip . we also obtain 
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d 



{^-l){B,-\{drV^{n)Y)- 



27 

1 



rdrMrs,^) (5.23) 



2j_ 
1 



Set Lw := a2V'o(^s) ■) for 02 defined in 

Since the constant C in (|4.1|) . (14.21) . (|5.14|) and (|5.21|) only depends on the data in sense of 
Remark 12.71 , combining ()5.22p with (I5.23|) , we conclude that the mapping L satisfies (jS.lOp 
with o{e) being independent of w. Thus we conclude D^Q = L. Since it is easier to verify 
()5.5p so we omit the proof. □ 

Corollary 5.3. Di,Q : C^^^Q/^-ji^) ~^ ^{~adA)^^^ ^ compact mapping. 

Proof. For any w G C^^^'^ ^^^(A), let ^/^o be a unique solution to (|5.19p . ()5.20p . Then there are 
constants C, C > depending only on the data to satisfy 

□ 

5.3. Local invertibility of V near C,q = {Id,ipQ ,Pq ,Vc)- For a local invertibility of V, by 
Lemma l5.H it remains to show that D^V is invertible and V is continuous in a Banach space. 

Note that for any a € (0, 1), ^^^^^(Co) and B^l{pc) are compact and convex subsets of := 
C2."/2(A/-,Mn) X C^'-/HK-+s) X C^'-^'iK^s) X C]:%^{A) and 'B(2) := ^^^i/laA) (A) re- 
spectively. We define, for z = (zi, Z2, z^, z^), 

I, II _n II , II |,(-a/2,9A) 

11^11(1):- l|2l||2W2X+F2|l3,„/2X-+||23|l2,a/2X-+,+ ll^4|li,„/2,A 

II II _ II ||(-"/2,9A) ^^-^^^ 
F4||(2) :- F4|li^„/2,A , 

and regard P as a mapping from a subset of *B(i) to a subset of *B(2)- 

Lemma 5.4. For any a G (0,1), V : B^\Co) ~^ ^Cu{Pc) is continuous in sense that if Cj 
converges to (00 in ^(i) o,s j tends to 00, then V{C,j) converges to 'P(Coo) in ^(2)- 

Proof. Take a sequence {Cj = (^j, fj)} in ^^^^■'(Co) so that ("j converges to Coo £ 

;Si^^(Co) in as j tends to 00. For each j G N, let Ej and £'00 be the solutions to the 
transport equation of ()4.8p with W* = W* and W* = respectively in ()4.9p where the initial 
conditions Ej{rs,'&), Eoo{rs,'&) in (I4.10p are determined by the data Cj and Coo respectively. 
Then, we have Q{Cj) = Ej\r^^ and Q(Coo) = -E^oolre^,- By the convergence of Cj in ^{1), 
W* converges to W* in c}'"^'^^-^ (7V+) n C^°^il p \(A/'+) as ? tends to 00. Because of 

{Cj} C Sct^^(Co)) by Proposition 13.11 and 14.11 {Ej} is bounded in C^^^j. ,)(-^rs) there is a 
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subsequence of {Ej}, which we simply denote as {Ej} itself, that converges to, say E^, in 
^{-a'/2r )(-^rt)- Then E'oo satisfies 

VF^ • DE^ = in M+, Eooirs,^) = Eo^{r„^) for all i9 G A, 

and thus E^o = Eoo holds in M^f. Moreover, this holds true for any subsequence of {Ej} that 

converges in C^l^^j2r )(-^rt)- Therefore, we conclude that Q{Cj) converges to Q(Coo) in ^(2) 
i.e., Q is continuous. One can similarly prove that TZ is continuous. By the definition of V in 
(j5.2p . we finally conclude that V is continuous. □ 

To prove the invertibility oi DyV : C^^^qj^^{A) — > C^f^^ gj^^{A), let us make the following 
observation. 

For a given w G Cj-f^gj^^iA), let V^""^ be the solution to (lETOD . ilQOl) . Define a mapping 

T by T : w I— )• — ^^/'^"'•'(rs, •), then T maps C^^^ q\^{^) into itself, and we can write DyV in 
Lemma 15.21 as 

D^Vw = ai{I -T)w. (5.25) 

By Corollary 15.31 T is compact. Then, by the Fredholm alternative theorem, either D^Vw = 
has a nontrivial solution w or D^V is invertible. In the following lemma, we show that D^V 
satisfies the latter hence D^V is invertible. 

Lemma 5.5. For a G (0, 1), D^V : C^{^f^^ q^^{A) — C"^'^^ q{^){^) is invertible. 

Proof. The proof is divided into two steps. 

Step 1. Let ai and 02 be defined as in (j5.8p . We claim 02 < and < 0. 
By Lemma [2121 we can write 02 as 

^2 _ "'•'^0 I _ 

Then, by (|23iD . HOG]) and (f3lB . we obtain 

(n - (r,)(i^o - (g^y^o (r,))^) 
02 = 7 < u. 

rsiBo - ^\Dip^{rs)\)~dr{ipQ - 'Po){rs) 



Since ai < is obvious from (|5.7p . there holds — < 0. 
Step 2. Claim. D^Vw = if and only if w = 0. 
If u; = then it is obviously D^'P'w=0. 

According to Remark |3.10| in spherical coordinates, (j5.19p and (j5.20p are expressed as 

dr{ki{r)dr'il^} + k2{r)Ag,i-iip = in Mj^,, drip — fJ-oi^ = on So, 

a I n r a I ^ T (5-26) 

Ou^V = on r^^,^^, drip = , , , on Tex 

fci(ri) 

where ki and /c2 are defined as in (j3.52p . 

Suppose DyVw = for some w G C'^l'^ aA)(^)- Then, by (I5.25p . there holds 

w{'d) = -^tjj(^){r„'d) for all ^ e A. (5.27) 
ai 
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We denote '0^'"^ as tp for the simpHcity hereafter. Let be a unit normal vector field on 
dA C Then is perpendicular to so w satisfies 

du^w = -—du^^{rs, •) = on dA. 
ai 

Then w has an eigenfunction expansion 

w = Cjr]j in A with Cj = / wqj 

3 = 1 ■'^ 

where every r]j is an eigenfunction of the eigenvalue problem 

— A§n-ir/j = XjiTfj in A, dy^rjj = on dA (5.28) 

with ||l2(^) = 1. We note that all the eigenvalues of (|5.28|) are non-negative real numbers. 

We claim that (j5.27p implies Cj = for all j G N. Once this is verified, then we obtain 
w = Q. 

For each € N, let us set wn '■= Z^jLi and 'i/^Ar(r, := X^jLi fo'^ one 
variable functions qj satisfying 

(fcig^.)' - = in (r„ n), (5.29) 

q'jirs) - fioqjirs) = 0, q'^{n) = (5.30) 

for j = 1, • • • , A^. Then ipNiTs-, ') converges to ip{rg,-) in -L^(A) so that (|5.27p implies 

qj{rs) = -Cj for all j G N. 

Suppose that cj 7^ for some j G N. Without loss of generality, we assume cj > 0. By the 
sign of ^ considered in Step 1, we have q'j{rs) < 0,qj{rs) < and q'jiri) > 0. We will derive 
a contradiction. 

(i) If Xj = then there is a constant niQ so that kiq'^ = tuq for all r € [rs,ri]. Since ki 
is strictly positive in [r^jTi], we have sgn niQ = sgn (?^(ri) and sgn niQ = sgn q'jirg). 
Then sgn q'j{rs) = sgn q'j{ri) must hold. But this contradicts to the assumption of 
q'j{rs) < < g^(ri). 

(ii) Assume \j > 0. By (j5.29p . we have 



\jk2qjdr = I {kiq^Ydr = ki{ri)qj{ri) - ki{rs)qj{rs) > 0. 

J Ts 

Then qj{rs) < implies max[rs,ri] Qj > so there exists ti £ [rs,ri] satisfying qjih) = 0. 
Also, by the intermediate value theorem, there exists t2 G ['^s,?"!] satisfying q'j{t2) = 0. 

If ti = t2 then the uniqueness theorem of second-order linear ODEs implies qj = 
on [rs,ri]. But this is a contradiction. If ti 7^ t2, assume ti < t2 without loss of 
generality. Then, by the maximum principle, qj = on [^1,^2] so qj = on [rs,ri]. 
This is a contradiction as well. 

Thus every cj must be hence w = 0. By the Fredholm alternative theorem, we finally 
conclude that D^.V is invertible. □ 
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Proposition 5.6. For any a G (0, 1), there exists a positive constant ai{< (73) depending on 
the data in sense of Remark \2. 7| so that whenever < cr < ai, for a given {^,(p-,p-,p, 
with <^i,<^2,^3 in Theorem[^ satisfying < Q < a for I = 1,2,3, there exists a function 
vl^{'^,ip^,P-,Pex) £ ^(-adA)(^^ ^^'^^ there holds 

V{^,yp-,P-,v:^)=Pex. (5.31) 

Moreover such a v*^ satisfies 

\Kx - ^c||S7„"r^ < + ,2 + ?3). (5.32) 
Proof Following the definition (lOnD . define, for the constant C in gH), V* : B* := S^^^Co) x 

V^{'^,iP-,P-,Vex,Pex) ■■= V{'^,iP-,P-,Vex) - Pex- 

Then, by Lemma [5.41 and the definition of above, for {(C~j)Pj)} C B* , if limj_j.oo |lCj ~ 
Cooll(i) + \\Pj -Poo||{2) = for some iCoo,Pinfty) G B* where || • |lfc=i,2 are defined in (|5.24p . 
then limj_j.oo ||P*(Ci)Pj) ~ 'P*{Coo,Poc)\\ = holds, i.e., is continuous in the sense similar to 
Lemma 15.41 Besides, by Lemma [221 'P* is Frechet differentiable at {Co,Pc)- In particular, we 
have D^V:t{Co,Pc) = DyV so it is invertible. So if a > is sufficiently small depending on the 
data and a then Lemma 15.11 implies (I5.3ip . (I5.32p is a direct result from Remark IA.21 □ 

Proof of TheoremUl Fix a G (0,1). By Proposition 15. 6[ for a given ^,92-, p^,Pex with 
?i5?2,?3 in Theorem [U satisfying < q < fii for / = 1,2,3, there exists a v*^ € C^^^ gj^-^{A) 
so that there holds (j5.3ip . For such a f*^. Proposition 13.11 and 14.11 imply that there exists a 
unique {p*,ip*,p*) in satisfying (I2.44p . (13. ip . Proposition [3Tl Lemma [42] and Proposition 
14.11 So {p, (p,p) = (p* ,P*)xj\f+ + (P-i ^-■•P-)Xj\f\M+ is a transonic shock solution satisfying 
Theorem ([T]) (a)-(d). Particularly, Theorem [T] (c) is a result from ()3.2p . ()4.ip and (I5.32p . Here, 
Xq = 1 in ri, and in Qp . □ 

6. Uniqueness 

6.1. Proof of Theorem [21 Fix C,^, := (^',(y9_,p_) and p^x with satisfying Theorem [1] (i)-(iii), 
and suppose that vi,V2 € C'(^1^9a)(^) satisfy 

V{C*,Vi)=Pex=V{(:*,V2). (6.1) 

If \\vj — fell 1 ak^^ ^ cTi in Theorem [H then by (j6.ip and Lemma [5.51 we have 

vi-V2 = -D,V-\V{C.,vi) - ViC,,V2) - D„V{vi - V2)). (6.2) 
If V were continuously Frechet differentiable, ()6.2p would imply 

II ||(-a,9A) / ^ II ||(-a,SA) 

SO that if o"! is sufficiently small, then vi = V2 holds, which provides Theorem [2] by Propo- 
sition 13.11 and 14.11 But it is unclear whether V is differentiable at other points near Qq = 
{Id, ifQ ,Pq ,Vc) for the following reason. 

For Vj{j = 1,2) chosen above, let W* be defined as in ()4.9p associated with C* = (^', (/?-,p_) 

and Vj, and let Ej € C^^^ p )(-^rt) t>e the unique solution to (j4.8p with (j4.10p . Then, Ei — E2 

satisfies M^^* ■ D{Ei - E2) = {W^ -W^)- D{E2 - E+) in with D = {dr, d^),, ■ ■ ■ , J by 
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Lemma 14.41 Then, by the method of characteristics used in section H?T| for any € A, Ei — E2 
has an expression of 

{El- E2)in,i}) 

r2ri—rs 

= [El - i^2)|x(2.,-r.;.,,^) -/ {W*2-Wl) •i^(i?2-i?o+)lx(t;n,^)dt (6-3) 

J ri 

' V ' 

where X{t;r,^) solves ([TOIl with W = W^. 

For Q in (j5.3p to be Frechet differentiable in sense of Lemma [5.2l we need to obtain a uniform 

bound of ^'"('-a 3A) foi' ^ sufficiently close to V2 with vi ^ V2-, but we only know, by 

ll''l-^2|li,,,A 

Lemmalia D{E - E^) G C°_„,r„)(-^rt)- 

But still, by (|4.15p and (|4.16p . it is likely that we have 

11-7- ||(l-a,9A) . ||(l-a,9A) 
< C\\vi-V2\\\^ 

for a constant C depending on the data in sense of Remark 12.71 Then (j6.2p may imply 

Ibi ~ '^alla a" '^'^^ — ^^ilbi ~ ^sIIq a"'^^^' ^'^'^ ^° ^1 ~ ^2 hold by reducing cji > in Theorem 
[T]if necessary. For that reason, we will prove the following lemma in this section. 

Lemma 6.1. For a G (^, 1), D^V defined in Lemma \5.S\ is a bounded linear mapping from 
^(i-a9A)(^) ^'^ itself and that D^V is invertible in C"^_^ q^-^{A) . 

Lemma 6.2. Fix a G (^, 1). Fix = (^f, (p^,p^) and vi,V2 G C^!^^ 9A)(^) fti '^2) defined 

in TheoremUl and := \\vj — Vc\\[ ^f^"^ ■ Then there is a constant C,(t^ > depending on 

the data in sense of Remark 2.1 so that whenever < ?2j S^4^\ ?4^^ ^ f^re satisfied there 
holds 

\\V{C.,vi) - V{C.,V2) - D,V{vi - V2)tX^^''^ 

^ ni I , V^ll ii(-«.9A)mi ||(l-a,aA) (^•4) 

< +?2 + 2^lbi -^^c|li,„,A )ll"l -^2||^,A • 

i=i 

Before we prove Lemma l6. II and Lemma 16.21 we assume that those lemmas hold true, and 
prove Theorem [2] first. 

Proof of TheoremlM Fix a G 1) and also fix ^,99_,p_,pez with ?i,?2,?3 in Theorem [T] 
satisfying < q < cr for a constant a G (0, fii] where ai is as in Theorem [TJ 

Suppose that there are two functions vi,V2 G C^!^^qa)(^) satisfying (15.3ip . (j5.32p with 
vl^ = Vj for j = 1,2. Then, by (|5.32p . (|6.2p . Lemma IHTTl LemmaE21 there holds 

11^1 -^2ii!:r''''^ < +,2 +^3)ibi -^2iii',r''^^ (6.5) 

Since C depends on the data in sense of Remark 12.71 there exists a constant a2 G (0,cji] with 
the same depndence as C so that we have 

C(ft + ^2 + <?3) < 3Ca2 < I. (6.6) 

Then, whenever < a < a2, from ()6.5p and (j6.6p . we conclude vi = V2- □ 

The rest of section [6] is devoted to the proof of Lemma 16.11 and Lemma 16.21 
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6.2. Proof of Lemma 16. IL Lemma l6.ll follows easily from: 

Lemma 6.3. Fix a € For F = (Fj)]^, G C^^_^^gg^^^Q^^JM+,m''), and g^g^ G 



c; 



(i-a,aA)(^)' 92 G C'fi-«,OSouar,,)(r^«,rJ; ^/^e /inear boundary problem 



dj{a^k{x,0)dki;) = djFj in M+ (6.7) 

Dtp -r - hqi/j = gi on Sq (6.8) 

{ajkix,0)dkip) ■ = g2 on T^^^, (6.9) 

{ajkix,0)dkip) ■ r = gs on Tex (6.10) 



/las a unique weak solution ip € C"'^'"(A/rt \ {9Sq U Srex)) H C"(A/r^). Moreover, there is a 
constant C depending only the data in sense of Remark \2. 7| so that there holds 

^ 11^ ||(1-",9A) , 11 ||(l-a,95oU9re,)^ (6.11) 

J=l,3 

The proof of Lemma 16.31 is given in Appendix [Bl 
Proof of Lemma \6.1\ By (|5.7p , for D^V to be well-defined in C^" (A) , the linear bound- 
ary problem (jS.lOp . (|5.20p must have a unique solution •0^"'^ G C2(AA+) n C°(A/'rt) for w G 
C(^_Q, g^)(A). This condition is satisfies by Lemma [6.31 Furthermore, (I6.1ip implies that D^V 
maps C(i_Q,aA)(^) itself, and it is bounded in C("_Q,aA)(^) ^ well. Also, the mapping 
w !->■ a27^(Co)V'^"'H^si ) is compact where 7^, 02 are defined in (j5.3p and (|5.8p . 

Suppose DjjVw^, = for some G C(i_q, aA)(^)' then by (j5.7p . we have 

7^(Co)a2 



It; 



Q(Co)«i 



V'^")(r„-)- (6.12) 



By (j6.1ip . ip^"^^ is in C"(A/'rt), then (j6.12p implies to* G C"(A), and again this provides = 
— 'P^^*\''"s, •) S C7"'^'"(A). This can be checked by arguments similar to the proof of 

Lemma 13. 5[ So one can follow the proof of Lemma 15.51 to show tu* = 0. By the Fredholm 
alternative theorem, we conclude that D^V is invertible in C"i_a dA)(^)- '— ' 

6.3. Proof of Lemma 16.21 To prove Lemma [6.21 we need the following corollaries of Lemma 
6.31 In this section, we fix ^, ip^,p^ with ?i, ?2 defined in Theorem [1] satisfying < q < a{l = 
1, 2) for a constant a > 0. 

Corollary 6.4. Fix a G (^,1). For j = 1,2, let (pj be the solution stated in Proposition \3.1\ 
with Vex = Vj satisfying q^"^ = \\vj — ■WcHiQ^yf^^ — ^ f^'^ ^ ^ (Oj^s]; '^'^^ with the transonic 
shock r = fj{x'). Set <j)j := {ipj — Lp^) o for G^. defined in (|4.3p . Then, there exist 
constant C,a^ > depending on the data in sense of Remark\2. 7| so that whenever {) < a < 



if < ?i, ?2 5 ?4"'^(j = 1,2) < fj are satisfied, then there hold 



ll-^i - '^211;:;^?"''^^^ < CH - V2tp'''\ (6.13) 
ll/i-/2llUr^<C^lh-^'2|&"''^^ (6.14) 
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Proof. Step 1. For defined in ([13]), let us write G-^{t,x') = {oj^'f\t, x'), x') then we have 

= t for all {t,x') G [r^n] x A. 
For convenience, we denote as Vj,u}j for v^^^\uj^^^\j = 1,2), then there holds 

vi{uji{t,x'),x') - V2{oJ2{t,x'),x') =Q for ah (t, x') G [r^, ri] x A (6.15) 
By the definition of Vj in ()4.3p . a direct computation shows that ()6.15p is equivalent to 

a;')]^i[Jf(i^^,) = - (t)2){t,x') + k[(l)j{t,x')x{ujj{t,x))]]=2 

where we set 

q{r,x') := k{ipQ(r) - (/9([(r) + V'_(r, x'))(l - x{r)) + rx{r) 
V'-(r,x') := ip-^{r,x') - ip^ir). 

By the choice of k in (|4.4|) . if (,2 in Theorem 12.41 is sufficiently small, then there is a positive 
constant go satisfying drq > (?o > for all r G [r^, ri] so that we can express uJi — 0J2 as 

, -/c(</.i-</>2)(i,x') + A:[</)j(i,x')xK-(t,x'))]}=2 

(wi - tJ2)(i,a; ) = -I — (6.16) 

Jq drq{aoJi{t,x') + (1 — a)co'2(t, x'))(ia 

for {t,x') G [rs,ri] x A. qq depends only on the data in sense of Remark 12.71 Furthermore, 
denoting as Gj, G"^ for G^^,G-j{j = 1, 2), since we have D{GjoGj^) = DGJ^ o DGjiGj^) = 
In for j = 1, 2, we obtain, in A/j^, 

[DG,{GJ^)] = -DGi{G^^)o{DG^^ -DG2^)oDG2{G2r^ . (6.17) 

with setting [DGj{G^'^)] := DGi{G^^)-DG2{G^^). 

Let us set Fens := dSo U dVex- Combining (j6.17p with (|6.16l) . we obtain 

li[^^^-(G7')]iiixr"^^ ^ ^(ii^i-'^2iiS;;^r)+ii^2(cui-a;2)ii;;;;Jr^ (6.19) 

(7) (7) (7) ( "1 (7) 

Step 2. For convenience, let us denote as aj^^ ,Ffr ,bi , fi^^' , gm for aki,Fk,bi, fif, Qm defined 
in dijD, (ISZD, (!3J5]1 and (l3:28D associated with ^,ip_,p_,ipj - f^Jj- Then, by dOD- dST]) . 
dSS]), and (|33i]l . for each j = 1,2, 0^ satisfies 

9fc(afcKy,O)9,0,) = a,.F^^Vafc([afcKy,0)-4?]5i<Ai)=:5fei?? in Ar+, 

9,0, - /io0, = #^ + (?^-^S'V^'/'i-(/io-/i-^^^'))0, =: on 

{dki{y,0)di())j) ■ = gi +{[aki{y,^)-aXi]di<l)j)-Vw =■ 02 ^w,rs 

{dM{y,0)di(t>j) -r = [vj - Vc) + {[dM{y,0)-d'-^i]di(t>j)-r =: g^^^ on Tex 

where dki{y,0), p^o are obtained from (I6.7p - (l6.10p through the change of variables (r, x') 
G^+{r,x'), and where d^^^ ,V"p , gm for A;,/ G {!,••• ttt, G {1,2} are obtained 

from a^j^i , bi\ Flf\ gli^ through the change of variables (r, x') i-7> Gj{r,x'), and thus they 
smoothly depend on D^, (p_,p_,99_), D(f)j and DGj{Gj^). 
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By ^B), IKlEh . (|33i]) . and by the definitions of d): ,Qm in (lOOl) . one can directly check 



I (1) (2)||(l-a,(95oU(9re^) ^ ^ e/ , , N 
|02 -32 lla,r„;, < CKj{(f>i,<P2) 

I (1) (2)||(i-Q,aA) ^,|| „ I ^ uc^ \ 

|03 - 03 IU,A < t^Ul^l - ^^2110, A + K*)d[(Pl,4>2) 



(6.21) 



where we set 



2 



_L /^.^ _L 



/ \ II / / ||(l-a,95oU(9rea;) 

Moreover, (|6.18|) and (j6.19p imply 

5{4>,A2) < ciu, - </'2ii;;;_^?"'"-^ + - ^2)<A2ii;;;_;^?"'"-^ 

By Proposition 13.11 and (|4.3p , each satisfies 



.22) 



ii^^-iiwS'""''"^^^ ^ ^(^1 + ^2 + 11^^- - -'^wtT^) (6-23) 



for a constant C depending on the data in sense of Remark 12.71 By (|6.18p . (|6.22p and (j6.23p . 

there exists a constant € (0,173] so that whenever < ?i,?2;?P'' < c'' are satisfied, there 
holds 

S{cl>ucP2) < - 02||S;"_^?'^'^^^^ (6.24) 

Here, the constants G and a'^ depend only on the data in sense of Remark 12. 71 Then, applying 
Lemma E3] to 0i - (/'2, by (fCT]) . ([6:21]) and (IHIMll . we obtain (f03]) reducing ct'' if necessary. 
Step 3. By arguing similarly as (|2.47p . and by the definition of <j)j, setting ip^ = ip^ — (p^ , we 
have 

U,-(r,,i?) - V-(/,-(i9),i?)]--2 
(/i - = ' ' ' ' (6.25) 

/; 9,((^o- - W + (1 - «)/2 

for X G A. By ([OS]) and dOSjl . we obtain 

ll/l-/2|ltaT<ai^^l-'A2||^;:^:^J"'"-^ + (?2+tf^ (6.26) 

We reduce cr^ further to obtain by ([OT]) and <^^. □ 



For any given a G (i, 1), let ip G C2(AA+) n C"(A/;t) be be a unique solution to (j6.20p with 
= 01 = 02 = and = vi - V2. 

Corollary 6.5. Fix a G (^,1). Let (f>i and 4>2 be as in Corollary \6.4i Then, there are 
constants C,a^ > depending only on the data in sense of Remark \2. 7| so that whenever 
< ?i, ?2 5 ?4^\ ?4^'* < o''^ O'fe satisfied, then there holds, 

II , , , |,(-o,9SoU9re^) , ^/ , , (1) , (2)m| ||(l-a,aA) fa n'7\ 
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Proof. By (I6.20p . and by the definition of ip, ip* := 0i — (/)2 — V' satisfies So (pi — (1)2 — u satisfies 

drip* - fiQip* = ^ - ^ on So 

{aki{y, 0)diip*) ■ = g^^^ - g^^^ on T^^^^ 

{aki{y, 0)diip*) ■ r = g^^^ - gf ^ - {vi - V2) on Tex- 

By the definition of g^^^ in (j6.20p . we emphasize that g3 := gg^'' — gg^^ — {vi — V2) satisfies 

llg.||(i-..aA) ^ ^^^^ ^ ^ ^(1) ^ ,(2)^^^^^^ 

where 5{(pi, <p2) is defined in (j6.22p . Then, by Lemma [Oj Corollary 16 .41 and (j6.2ip . we conclude 
that if is chosen sufficiently small depending on the data in sense of Remark 12. 7| then 
tp* = (pi — (p2 — Ip satisfies 

II , n{-a,dSoVJdV,,:) ^, 1 ^(1) _L^(2)^|^, „ |[(l-"-9A) 

ll^*"l,«,A/'+ < + ^2 + +?4 )\\Vl-V2\\a^A {b.2,'6) 

for a constant C depending only on the data in sense of Remark 12.71 □ 
Fix in 

"fo ,Po )• For simplicity, let us denote as V{v) for ip^,p-,v) 

Proof of Lemma \ 6.2[ Step l.For convenience, for each j = 1,2, let us denote as Vj, Qj and 
TZj for V{C*,Vj), Q{C*,Vj) and Tl{C*,Vj) respectively where Q,TZ are defined in (|5.3p . Then we 
can write 

V1-V2- D.,V{vi - V2) = hTZi + I2Q2 + h (6.29) 

with 

ii = Qi-Q2- D^Q{vi - V2), /2 = 7^l - 7^2 - D^n{vi - V2), 

h = (Q2 - Qo)D^n{vi - V2) + [Ui - no)D„Q{vi - V2), 

where we denote Qq = Q{Id,ipQ ,Pq ,Vc), TZq = Tl{Id,ipQ ,Pq ,Vc)- We estimate liR.i-,l2Q2-,h 
separately. 

By Lemma l6.H Proposition 13. 11 (I4.14p and (I4.15p . we obtain 

l|, ||(l-a,aA) I ^(1) _L^(2)\|L, „ ||{l-°.aA) I . 



Step 2. Let ip G C2(7V+) n C°(A/;t) be a unique solution to ([5T9]) . (fOO]) with w = vi-V2, 
then by (|6.3p and the definition of Q in (j5.6p , we have 

ll{'d) = {Ei-E2 -a2ip)\x{2ri^s-,Ti,i))+a2{ip\x{2ri^s-,Ti,i))-'<P{rs,'&))+I*{^) (6.31) 

where 02,/* are defined in (|5.8p . (|6.3p respectively. First of all, by (j4.14p . (|4.15p . Corollarv 16.41 
and (|6.18p . we obtain 

II r ||(l-a,9A) . , , (2)x|| ||(l-a,9A) Qn\ 

P*||'a <C{(.i + (;2 + <^l')\\vi-V2\tf, '. (6.32) 
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Secondly, for Vq defined in Lemma 14.31 let us set Wq := so that we have Wq = 
(1, 0, ■ ■ ■ , 0) in (r, -d) coordinates which is a spherical coordinate system we specified in section 
m Then, by Lemma |4.3| there holds 

\\W; - \\[~";^+''^^'''^ < C{<; + + ?f ) (6.33) 

for a constant C depending on the data in sense of Remark 12.71 where Wi is defined in the 
paragraph of ()6.3p . Thus, by (|6.33p and Lemma [Oj we obtain 

||V'|xi{2ri-r-s;ri,-) " "0(^5, O lla,!"'^'^'' 

a, A 



< C(?i + c;-2 + Oll^i - ^2" 



for a constant C depending only on the data. 

To complete step 2, it remains to estimate ||(£'i — E2 — a2ip){rs, Olla a"'^"^^ because Xi(2ri — 
r^; ri, I?) lies on Sq = {r = r^} H M for all ■(9 G A. 

By the definition of 02 in (|5.8p . 02'?/' is expressed as 

a2ipirs, •) = a^^Vl^'s, •) + ■) (6.35) 

with 

-(i^o-i(9.^^(r.))^)^5.(v^o--^j[)(^.)' ^ ^ 
where 03 is defined in ()5.23p . So, we need to decompose Ei — E2 into two parts so that one is 
comparable to a^'^t/j, and the other is comparable to a'^^ip. 



with 



For each j = 1, 2, by (j4.10p . we have 

(Mi 

(So - i|Z)^-i(^)V(^j|2)^ ' U. 
V((^_ -(^j)|(Z)^-i)^(D*-i 



\D^f-^Viip^-ip,)\ 



where every quantity is evaluated at (/j (-!?), -i?). Here, ipj and fj are as in Corollary 16.41 and 
z^P'' is the unit normal on Sj = {r = /j(i9)} toward the subsonic domain A/j" of ipj for j = 1, 2. 

Set f/o(r) := {Bq — ^\dripQ {r)\'^)^ . Then, for each j = 1,2, we can write Ej{rs,'&) — 
E+ {fj (t?)) = ^ - ^ + ^P^.o ^here ps,o, ^o" is defined in (|233|) and LemmalHrespectively, 
and where all the quantities are evaluated at (/j(i?), i?). Then, by Lemma |4.4| we can express 
{Ei-E2){rs,^) as 

{Ei-E2)irs,^) = zi + Z2, (6.38) 

p^^ Ps,0{m) .l ^ _ r (Ps,0-P0+)(/.(^)) il ,q^ 
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Here, [Fj]]^^ is defined by [Fj]]^2 ■= - -^2- 

We claim that, for a constant C depending on tlie data in sense of Remark 12.71 there hold 

-a?V(r„.)||l^,r''^^ < C,A\vi-v,tX^'''\ (6.40) 
\\z, - a^^^in, Olll^^r-^^) < C,4v, - v.tX^'^'K (6.41) 

where we set := ?i + ?2 + S'l"'^^ + '^f^- 

Step 3. (Verification of ([OT]) ) We rewrite Z2 in (f09]) as 

_ [(Ps.0-P^)|r=/,W]]=2 . 1 1 Nil 

By ()2.33p . Proposition 13.11 and Corollary ()6.4p . we obtain 

II (2)||(l-a,aA) ^ II ||(l~a,9A) . .^ 

By (j6.25p . we can write 

z«(^) - a«^(r„^) = (^) (^s , ^) - ^- (/, (^) , nU - ^i^s , ^) 
^ ^ Uo{rs) 

with 

-(1),^. ^ /o dr{ps,o-p^){aM^) + (1 - a)/2(^))rfQ 

C/o(r,) J^driifo - ^t){ahm + (1 - a)/2(^?))(ia' 

where each cpj is as in Corollary 16.41 and where we denote as "0- = V- ~ V'o • By Corollary 16.41 
and Corollary 16. 5^ there holds 

The constants C in (j6.43p . ()6.44p depend only on the data in sense of Remark 12. 7i Then, (j6.4ip 
is obtained by (I03D and dOIl) . 

Step 4. (Verification of iK^ ) By and (IS22D, it is easy to see that for 4^V(?'s, •) 

in (|6.35p . there holds 

4'V(rs, •) = lPUrs)dript{rs) q^^^,^^ .) = a,dri^{r,, ■). (6.45) 

(^_l)(i?o-i(5.v.+ (r,))2) — 

Similarly to step 3 aboye, we rewrite zi in ()6.39p as 

^, -„+^^ul 1 11 I r.pi'^ ^ I Po(^s)-Ps,o(/,(^)) .i 

'""^"^^^^C/; C/o(/,(^))^^=' + ^^C/, Uf^ UoUm (6.46) 



where we set 



U*:= (Bo-\\V^ifjm,m^'- 
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By (IPs]) . Proposition EH (fOTD and Corollary [631 zf^ in ([06]) satisfies 

II (2)||(l-a,9A) , ^ II ||{l-a,9A) /lv^ 

Pi^ll^ ' < C?*||t;i -z;2||„,A ^ (6.47) 
For each j = 1,2, we can write 
+ r 1 1 ^ 



7P^(^.)/?l^H^)-V(¥^i-V^o)l(/,W,tf) (6.48) 

with f3i'\^) = drift if + iV(^, - (/^o^)|(/,W,^). 

By Proposition 13. 1[ we have, for a constant C depending on the data in sense of Remark [2 .71 

\\a?-a,f\\[;^^f^Uc,^. (6.49) 

So, if we write z^^ — a^dri^irs, •) as 
zj^^ - a4drip{rs,-d) 

= a4([9r(¥'j-V?o )l(/,(tf,^))]]=2-V'(^^s,i?)) + [(aJ'^-a4)-V((/?j-v9fj")|(/.(^),i»)]]=2, 

then, by Proposition 13.11 ([4.3p . Corollary 16.41 Corollary 16.51 ([6.49p . there is a constant C 
depending on the data in sense of Remark 12.7! so that there holds 

- a,dMrs,mtA"''''^ ^ - v^tp'^K (6.50) 

Thus ([6.40p holds true. Finally, ([6.40p and ([6.4ip imply, for a constant C depending on the 
data in sense of Remark I2.7[ 

Similarly, one can also show 

\\Q2l2tp'''^ < C,.\\vi - V2tp'^^. (6.52) 
As (I6.52P can be checked more simply than (I6.5ip . we omit the details for the yerification of 

Finally, combining (16.30p . (16.5ip and (I6.52p all together, we obtain (16.4p . and this competes 
the proof of Lemma 16.21 

□ 

Appendix A. Proof of Proposition 15.11 

Before we prove Proposition 15. H we first prove the following lemma. 

Lemma A.l (Right Inverse Mapping Theorem). Let be a compactly imbedded subspace of 
a Banach space Also, suppose that €2 is a subspace of another Banach space '^2- For a 
given point (xo,yo) G G^i x <t2, suppose that a mapping % satisfies the following conditions: 
(i) 1-i maps a neighborhood of xq in 5Si to 5^2, and maps a neighborhood of xq in Ci to 
€2, with 7i{xo) = yo, 
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(ii) whenever a sequence {xk} C near xq converges to x* in the sequence {Hixk)} C 
£2 converges to T-L{x^) in 

(iii) T-LjUS a mapping from *Bi to also as a mapping from £1 to £2? is Frechet differen- 
tiable at xq, 

(iv) the Frechet derivative DT-LIxq), as a mapping from QSi to ^2, o-lso as a mapping from 

to (^2, is invertihle. 

Then there is a small neighborhood V(xo) of xq in Ci, and a small neighborhood yV(yo) of 
in C2 so that % : V(xo) W(?/o) has a right inverse i.e., there is a mapping 'H~j^gy^^ : W(yo) — ^ 
V(xo) satisfying 

^-^m,ht = id (A.l) 
and moreover T-ij^gj^^ satisfies T-L'^^^^ o 7i{xo) = xq. 
Proof. Denote D'H{xq) as DT-Lq- For t € M+, and v,w G Ci, define 

n*{t, V, w) := ^ [n{xQ + t{v + w)) - n{xo)) - DUqw. 

For r > and G £1, set Br{ao) := {y G Ci : \\y — aollci — '^l- ^'^^ ^ fixed (t, € 
(0,eo] X (0), we define 

niJv) ■.= Dno\n*){t,v,w)-v. 

If the constant Eq is sufficiently small then we have Tily^ maps -Bi(O) into itself. Since -Bi(O) 
is bounded in Ci, -Bi(O) is a compact convex subset of By Lemma [A.ll (ii) and (iii), 
the Schauder fixed point theorem applies to Ht^^. Thus ^j^, has a fixed point v^^ = V:^{t,w) 
satisfying 

By the definition of 7^*^^,, f* satisfies 

^n{xo + t{v^ + w)) = n{xo) + tDnow. ^ 

We claim 

for z = 'H{xq) + tWHovu. Obviously, T~iJ^^gj^^{z) in ()A.3p satisfies (jA.ip . So it remains to check 
if T-L]ilght in (|A.3p is well-defined. 

Set W(yo) := {yo + tDHow : t G [0,eo],''i' G 5Bi(0)}. Since DUq is invertible, W(yo) covers 
a small neighborhood of yo in ^2- Suppose yo+tiDT-LQWi = for ti, t2 G [0, Eq] and 

tt?2 G dBi{0). If = then ^2 must be as well. If 7^ then we have DT-LqWi = ^DT-LqW2 
and this implies vui = ^W2. So we obtain ti = t2 and wi = W2. 

By (|A.3p . we also have 7^^]^;^^ o •^(xo) = HRlghtC^i^o) + Ovu) = xq. □ 

Remark A.2. By (jX3]) . i/i ere /loZds 

p-H^) - xoU, <2t = 2^— < 2||Z??^o-i|l||z - 'H{xo)\W, 

where we set WDT-Lq^W := sup||^||^_^=o il-^'^o ^Iki- 
Proposition [STT] can be easily proven from Lemma lA.li 
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Proof of Proposition \5. il Let us define T-L{x,y) by 

n[x,y):={nx,y),y)- (A.4) 

In LemmaEIl replacing «Bi by *Bi x «B2 , by £1 x £2, ®2 by QSs, £2 by £3, xq by {xq, yo) 
and yo by (0, ?/o)) ^ in (|A.4p satisfies Lemma lA.ll (i)-(iii). In particular, the Frechet derivative 
of 7i at {xo,yo), which we write as DHq, is given by 

DHoix, y) = {DxTqx + DyToy, y) 

where we denote -D^^, J-'(xo, yo) as (D^J^o, DyJ^o). By Proposition 15. II (ivl. Ti satisfies Lemma 
lA.ll (iv). Therefore, by Lemma lA.ll T-L has a right inverse 'H]^lgj^^ in a small neighborhood of 
{xo,yo) in '^^i 2^2- In particular, for y € £2 with ||y — yo||c2 being sufficiently small, 'H'j^gj^^ is 
well defined. Write T^Righti^^y) — ix*iy)jy) then we obtain T{x*{y),y) = 0. □ 

Appendix B. Proof of Lemma 16.31 

Proof of Lemma \6.3[ 

Step 1. Fix F* = {F,)-^, G C^ri-a,asouar..)(-^^ ' I^")' ^^^^ ^LSs* e C("^_,^(A), 52* e 
Ca-a,a5ou9re.)(r»,rJ- Then, there are sequences {F^^) = {F^^^)"^^^} C C"(AA+,M"), and 
rf^lJfff } C C7-(A), {^f ^} C ^"(1^) satisfying 

lim ||^(fe)_^*||(i-",a5ouar..)^ ^.^^ II^W_^*||(i-«,aA) 

fc— )-00 (^y-Nrs fc— >00 ' 



lim II o^'^) - .*||(l-".95oU9re.) _ II (fc) _ *||{l-a,9A) _ . 

nm iig'a 52 110,?^^^ - nm 11^3 5'3lla,A - ^• 

fc— fOO ' ' " fe— >-CXD ' 



For each k, let V^''^ S C2(7V+) n C°(A/;t) be the unique solution to (16771) - (16301) with F = 

F^^\gi = g\^^ for / = 1,2,3. The unique existence of il^^^^ can be proven by the method of 
continuity as in the proof of Lemma [ 



If (I6TT]) holds true for F e C"(A/;t, M"), 5i, 52 e '^°(A),52 S C°(r^,^J, then we have 

lim IIV^C^) -V'(-)||(-"f^^°^^^=-) =0. (B.l) 

fc,m— >00 liC'iA/r-s 

By (|BT]) . V^'') converges to a function ^* in for any compact subset of M^^, so -f/;* satisfies 
()6.7p with F = F* in 7V"+. Also, V* satisfies ()6.8p - ()6.10p on the relative interior of each 
corresponding boundary with gi = for / = 1,2,3. Since {ip^'^^} is a Cauchy sequence in 
Cl-a,ds,udr...)(^^s), we have r e Cl:^a,9s,udr.^)i^^J ^ C-{M)- Then, by |LlIl Corollary 
2.5], iiQ > in dSSD implies that V* G C'2(AA+) n C°(Ar;I"J is a unique solution to (fHTTP - dUTUIl . 
Moreover, ij)* satisfies (16. lip by (jB.ip . Thus, it suffices to prove (I6.11j) for 



F G C°(A4t,K"), 51,52 G C"(A), 52 e C"(r^,,J. (B.2) 

So, for the rest of proof, we assume (|B.2p . 

Step 2. By the local scalings as in the proof of Lemma 13. 5( one can easily check that if 
G C2(A/'+) n C°(7Vrt) solves (f6Tl) - (|6.10p . then there is a constant C depending on AA^ and 

a so that there holds 

,{-a^5ouar..) ^ c-dV^U^^^ + if(V',F,5i, 52,53)) (B.3) 
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where we set 

T^/ I T? \ II p||(l-a,a5oU9re^) (l-a,95oU9re^) , n ||(l-a,aA) 

/=1,3 

So it suffices to estimate IV'Iq ^ Ar+ • 

Since A/"^ is a cylindrical domain with the cross-section A in (r, x')— coordinates in (I2.22p . 
there is a constant > depending only on n, A to satisfy, for any xq € A/"^"^ and i? > 0, 

1 vol{BR{xo)r\M+) ^ 

— < ^ < Kq. 

Kq it" 

So, if there is a constant i?^, > and M > satisfying 

pn-2+2a / jl^V' l' < for ah i? G (0, i?,] , (B.4) 

then there holds 

-"'* 

assuming i?* < 1 without loss of generality where the constant C depends only on n, a, A, r^, ri. 
So, we devote the rest of proof to find a constant M > satisfying ()B.4p . 

To obtain ()B.4p . we need to consider the three cases: (i) Bji{xo) C M,^^, (ii) xq € 9A/^+ \ 
{dSo U dTex) and Bji{xq) n (SS'o U SFex) = 0, (hi) G SS'o U c^Fe^. More general cases can be 
treated to these three cases. Also, case (i), (ii) are easier to handle than case (iii). So we only 
consider case (iii). 

Step 3. Assume xq G dSo because the case of xq G dT^x can be treated more simply. For a 
fixed point xq G OSq, let xq = (rs,a;g) in the (r, x')-coordinates given in (j2.22p . Then, there is 
a constant Rq > depending on A, and a smooth diffeomorphism h defined in a neighborhood 
Or^'^ of Xq in A so that h flattens dA near x'q, and moreover the followings hold: 

(a) h{x'o) = G (b) for any x' G n ^A, h{x') G M""^ x {0}, 

(c) for any (r,x') G 5r„(xo) nA^^, (r,/i(x')) G (r^n) x M"-^ x M+. 
For a constant R > and yo := (^s) ^(^o)) = (^si 0) G M", let us set 

:= Bnivo) n {(r,/i(x')) : {r,x') G B^j^ (xq) n A^^}, 
Sr„,fi:=Si?(yo)n{(r,/i(x')) : {r,x') e Br,{xo) nVy,}, (B.6) 
S5o,i?:=^i?(yo)n{(r,/i(x')) : (r, x') G ^^(.(xo) n Sq}. 

For (r, x') G A/"^"!^, let us write y = {r,y') = (r, /i(x')), and (p{y) = ip{r,x'). By (|6.7p - (|5.9p . (j) 
satisfies 

dkiaki{y,0)di(l)) = dkFk in Pr, dr(j) - Ho(j) = gi on Ssq^r 

(B.7j 

{aki{y,0)di(t>) -Uw =92 on Sr^,/? 

where aki, Fk, gi, g2 are obtained from aki, Fk, gi, g2 in (j6.7p - (|6.10p through the change of vari- 
ables (r, x') (r, h{x')). So ^=1 is strictly positive, and the regularity of aki, Fk, gi, g2 are 
same as the regularity of aki, Fk, gi, g2 in T^R- Here, is the inward unit normal to T^r^^R- 
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Since h is smooth, and A is compact, there is a constant i?i > depending on n, A so that we 
have 

dV2R, \ (S5o,2i?i U Sr„,2Rj C dB2RAyo)- 
step 4. Fix R G (0, Ri], and set ^dB,R ■= dVn n dBniyo). Write (j) in ([RT]) as (/> = u + 
for a weak solution u to 

9fc(aM(yo,0)9iu) = in P/j, u = ()) on T^qb^r 

drU - fXQU = on Y.So,R, {akiivo, 0)diu) • = on Sr^,/?. 

Such u G uniquely exists by Lemma 12.51 and by a basic estimate for harmonic 

functions as in |HLl Lemma 3.10], there is a constant C depending on the data in sense of 
Remark 12.71 to satisfy 

/ \Du\'^dy < / \Du\'^dy for < qi < g2 < R. (B.9) 

Step 5. Then, by ([BTI) and (iRSl) . w = (f) - u satisfies 

dk{dki{yo,^)diw) = dk{[dki{yo,G)-aki{y,G)]di(t) + Fk) in Vr, 
drW - flow = gi on T.So,r, 

{akiiyo,0)diw) ■ = {[dkiiyo,0) - dkiiy,0)]di(p) ■ + g2 on Sr^,^, 
i(; = on Sas^R. 

By (j3.41|) . (j3.52|) and (jB.lOp . there is a constant A, C > depending on the data in sense of 
Remark 12.71 so that there holds 

X \Dwfdy<C{ \hi\\Dw\dy+ \h2\\w\ dAy+ \hs\\w\dAy) (B.ll) 
where we set 

\hiiy)\=6*iyo,y)\Dc^{y)\ + \F{y)\ 

\h2{y)\ = 6*{yo,y)\D(t>{y)\ + \F{y)\ + \w{y)\ + \h{y)\ (B.12) 

\h3{y)\ = 6*{yo,y)\D(t>{y)\ + \F{y)\ + \~g2{y)\ 
with 5*{yo,y) = (EL=i \m{yo,0) - dki{y,0)\^y/^ . 

Before we proceed further, we first note that, by the definition ol (j),F = {Fk)^^^, gi, g2 in 
(|B.7p . for any /3 G [0, a], there is a constant C(/3, A) > depending on (3, A so that we have 



(-«,Sso,-RnSr„,H) ^ r'fR AMU/,ii(~"''5'^ouare:r) 

I^ll5*'^r^''^^^^(/5,A)||F||(;;^'^^«^^^-) with F = (F,)Li 
|5i||/3,S5o,/?(^-"'^^^) < C(/3,A)|bi||gX"''^) 



(B.13) 



Let us set d{y) := dist{y, '^So,R ^ ^r^.i?) and 

M2{R) .- Mi[R) + \\gi\\Qj^ + ll52|lo,Sr„,jj 



(B.14) 
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Using (IB.12|1 . (|B.13p and the continuity of dki{y, 0) in 5*{yo, y), we get, for any y G Vr^^, 
\h,iy)\<M,iRMy)]-'+", 

\\h2{y)\ - \w{y)\\ < M2(i?)[d(2/)]-i+", \h3{y)\ < M2(i?)[%r ^+". 
By the Poincare inequahty with scahng, and ()B.15p . for any e > 0, we have 

\h2\\w\dAy 



(B.15) 



i 

<C{e,n,a)R[M2iR)f f [d{y)f^-^+'^UAy + {e + C{n)R)f \Dw\'^dy. 

•'So,, B •'Cr 



(B.16) 



By (|B.6p . for any R < Ri, S^q,/? satisfies 



^So,RC{{r,y[,--- , y^^i) : r = r„ y^^^ G {0, R), {y[, ■ ■ ■ .yU) ^4" '^0)1 (B-17) 



where B^{0) is a ball in M of radius R centered at G M . This provides 

f-R 

[d{y)]'^(-^+^)dAy <C{n)vol{B)!^~'^\0)) / f^-^+'^Ut<C{n,a)R''-^+^''. (B.18) 



Note that a > ^ is an important condition in (]B.18p . By (IB.16P and (IB.lSp . we obtain 
/ \h2\\w\dAy <Cie,n,a)R''~^+^''[M2iR)f + ie + C{n)R) f \Dw\'^dy. 

One can treat \hi\\Dw\dy, ^ \h3\\w\dAy similarly, and combine all the results to- 
gether so that (IB.lip implies 

\[ \Dw\^dy <C{e,n,a)R''''^^'^'^[M2{R)]^ + {e+C{n)R) f \Dw\^dy. (B.19) 

Choose e = and reduce Ri if necessary so that (jB.lOp implies 

/ l-Dwpdy < C(n, a, A)i?""^+^"[M2(i?)]^ for all i? G (0,i?i]. (B.20) 
Step 6. By (jR9]) and ^(K20\i . whenever < Q < R < Ri, (p satisfies 



/ m^dy < C((|)" / m^dy + i?"-2+2"[M2(i?)]^ 



\-iy(p\ uy -t- n [iv-i2i^n/^ 
So, by |HH Lemma 3.4] and (jB.lSp . we obtain 

?2-2a 



^ m'dy < CR--'+'^{^\m\l,^^^^ + [M2{R)f) 

< Ci?"-2+2"((||V'|i;-^f°^^^-^)2 + [M2{R,)f). 



(B.21) 



We note that the estimate ()B.2ip holds true for any choice of xq = h^^{yo) G and any 
R G (0, Ri], moreover Ri is uniform for all xq G M^^. Since A is compact, and dA is smooth, we 
can take a constant R* > depending on A so that, by ()B.2ip . the solution ip to ()6.7p - ()6.10p 
satisfies 

/ IZ^V'P <OT-2+2"(||V.||(-°'^f^^r..)^^^^2 (B_22) 
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With Ao - +I.z=i,3ll5i|lo,A +ll52|lo,r„,.. 

By the interpolation inequaUty in \GT\ Lemma 6.34], for any e > 0, we have 

So, assuming R* <\ without loss of generality, by (|B.3p . (jB.Sp . we obtain 

||^|j(-a^5ouar..) ^ +i^(F,5i,52,53)) 

where K{F, gi, g2, gs) is defined after (jB.Sp . 

Since /^o > in (|6.8p by Lemma [231 by the uniqueness of solution( |LIH Corollary 2.5]) for 
(|6.7p - (|6.10p . we have |^lo_^^+ ^ CK{F, gi, g2, ga). The proof is complete. □ 



Acknowledgments. The authors thank Marshall Slemrod for helpful comments on this 
work. Mikhail Feldman's research was supported in part by the National Science Foundation 
under Grants DMS-0800245 and DMS-0354729. 

References 

[AM] Amann, H., Ordinary differential equations: An introduction to nonlinear anaJysis(translated in 

English), Walter de Gruyter, (1990). 
[BR] Berard, Pierre H., Spectral Geometry:Direct and Inverse Problems, Springer- Verlag, (1986) 

[CKK] Canic, S., Keyfitz, B. L., and Kim, E. H., Free boundary problems for nonlinear wave systems: 

Mach stems for interacting shocks, SIAM J. Math. Anal. 37 (2006), 1947-1977. 
[CKL] Canic, S., B. Keyfitz, L. and Lieberman, G., A proof of existence of perturbed steady transonic 

shocks via a free boundary problem, Comm Pure Appl. Math., 53 (2000), 484-511 
[CGF] Chen, G.-Q., Chen, J. and Feldman, M., Transonic Shocks and Free Boundary Problems for the 

Full Euler Equations in Infinite Nozzles, J. Math. Pures Appl, 88 (2007), no.2, 191-218. 
[CFl] Chen, G.-Q. and Feldman, M., Multidimensional transonic shocks and free boundary problems 

for nonlinear equations of mixed type, J. Amer. Math. Soc. 16 (2003), 461-494. 
[CF2] Chen, Gui-Qiang; Feldman, Mikhail Steady transonic shocks and free boundary problems for 

the Euler equations in infinite cylinders. Comm. Pure Appl. Math. 57 (2004), no. 3, 310-356. 
[CF3] Chen, G.-Q. and Feldman, M., Existence and stability of multidimensional transonic flows 

through an infinite nozzle of arbitrary cross-sections. Arch. Ration. Mech. Anal 184(2007), no 

2, 185-242. 

[CHI] Chen, Shuxing Stability of transonic shock fronts in two-dimensional Euler systems. Trans. Amer. 

Math. Soc. 357 (2005), no. 1, 287-308 
[CH2] Chen, S.-X., Transonic shocks in 3-D compressible fiow passing a duct with a general section for 

Euler systems Trans. Amer. Math. Soc. 360(2008), no. 10, 5265-5289 
[CHS] Chen, S.-X., Compressible fiow and transonic shock in a diverging nozzle Commun. Math. Phys. 

289(2009) 75-106 

[CY] Chen, S.-X., and Yuan, H.-R., Transonic shocks in compressible fiow passing a duct for three 

dimensional Euler system. Arch. Ration. Mech. Anal. 187(2008), no 3, 523-556 

[CFR] Courant, R. and Friedrichs, K. O., Supersonic Flow and Shock Waves, Springer- Verlag: New 

York, 1948. 

[GT] Gilbarg, D. and Trudinger, N., Elliptic Partial Differential Equations of Second Order, 2nd Ed., 

Springer- Verlag: Berlin, 1983. 
[HL] Han, Q. and Lin, F., Elliptic partial differential equations, Courant Institute of Math. Sci., NYU, 

1997. 

[HS] Hestense, M.R. Calculus of Variations and Optimal Control Theory Wiley & Sons, New York- 

London-Sydney, 1966 

[KU] Kuz'min, A.G., Boundary Value Problems for Transonic Flow John Wiley & Sons, 2002 



54 



MYOUNGJEAN BAE AND MIKHAIL FELDMAN 



[LXY] Li, Jun; Xin, Zhouping; Yin, Huicheng On transonic shocks in a nozzle with variable end pres- 

sures. Comm. Math. Phys. 291 (2009), no. 1, 111-150. 

[LIl] Licbcrmau, Gary M., Local estimate for subsolutions and supcrsolutioiis of oblique derivative 

problems for general second order elliptic equations Trans. Amer. Math. Soc. 304(1987), 343-353 

[LI2] Lieberman, Gary M., Holder coninuity of the gradient of sotluions of uniformly parabolic equa- 

tions with conormal boundary conditions, Ann. Mat. Pura Appl. (4), 148(1987), 77 99 

[LI3] Lieberman, Gary M., Oblique derivative problems in Lipschitz domains: L Continuous boundary 

data Boll. Un. Mat. Ital., (7)1-B(1987), 1185-1210 

[LI4] Lieberman, Gary M., Oblique derivative problems in Lipschitz domains: ILDiscontinuous bound- 

ary data J. reine angew. Math., 389(1988), 1-21 

[LU] Liu, T.-P., Nonlinear stability and instability of tranosnic gas flow trhough a nozzle. Comm. 

Math. Phys., 83 (1982), 243-260 

[LY] Liu, L. and Yuan, H.R., Stability of cylindrical transonic shocks for two dimensional steady 

compressible Euler system, J. Hyperbolic PDE, 5(2008), no 2, 347-379 

[OO] Ockendon, H. and Ockendon, ,L R., Waves and Compressible Flow, New York : Springer, 2004 

[SZ] SZILAGYI, T., Inverse function theorem- without assuming continuous differentaibility 

Ann.Univ. Sci. Budapast Eotvos Sect. Math., 20 (1977), 107-110 

[XYl] Xin, Z.-P. and Yin, H.-C, Transonic Shock in a Nozzle L Two Dimensional Case, Comm. Pure 

Appl. Math, 58 (2005) ,999-1050 

[XY2] Xin, Z.-P. and Yin, H.-C, The transonic shock in a nozzle, 2-D and 3-D complete Euler systems, 

J. Differential Equations, 245, (2008), 1014-1085 

[YUl] Yuan, Hairong On transonic shocks in two-dimensional variable-area ducts for steady Euler 

system. SIAM J. Math. Anal, 38 (2006), no. 4, 1343-1370 

[YU2] Yuan, H.R., A remark on determination of transonic shocks in divergent nozzle for steady com- 

pressible Euler flows Nonlinear Anal. Real World Appl. 9 (2008), no 2, 316-325 

M. Bae, Department of Mathematics, Northwestern University, Evanston, IL 60208, USA 
E-mail address: bae@math.northwestern.edu 



M. Feldman, Department of Mathematics, University of Wisconsin, Madison, WI 53706-1388, 

USA 

E-mail address: feldman@math.wisc.edu 



